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Atom interferometry is now reaching sufficient precision to motivate laboratory tests of general 
relativity. We begin by explaining the non-relativistic calculation of the phase shift in an atom 
interferometer and deriving its range of validity. From this, we develop a method for calculating 
the phase shift in general relativity. Both the atoms and the light are treated relativistically and all 
coordinate dependencies are removed, thus revealing novel terms, cancellations, and new origins for 
previously calculated terms. This formalism is then used to find the relativistic effects in an atom 
interferometer in a weak gravitational field for application to laboratory tests of general relativity. 
The potentially testable relativistic effects include the non-linear three-graviton coupling, the gravity 
of kinetic energy, and the falling of light. We propose specific experiments, one currently under 
construction, to measure each of these effects. These experiments could provide a test of the principle 



■ metrics including the Lense-Thirring effect, the expansion of the universe, and preferred frame and 
location effects. 
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I. 



INTRODUCTION 



Atomic physics has made rapid progress recently, yielding several high precision techniques and experiments. Atomic 
clocks can remain synchronous to 16-decimals [1] and have been used to search for the time variation of the fine 
structure constant [2-4]. In particular, atom interferometry has become a high precision tool that has been used 
successfully for a variety of applications including gyroscopes [5], gradiometers [6, 7], and gravimeters [8, 9]. In light 
of this, it is interesting to consider the possibilities for using atom interferometry to test general relativity in the 
lab. Indeed, the first atom interferometry experiment to push the limits on the Principle of Equivalence (PoE) [10] is 
currently under construction. 

Experimental tests of general relativity (GR) have gone through two major phases. The original tests of the 
perihelion precession and light bending were followed by a golden era from 1960 until today (see e.g. [11]). These 
tests were in part motivated by alternatives to Einstein's theory, such as Brans-Dicke, designed to incorporate Mach's 
principle [12]. More recently, the cosmological constant problem suggests that our understanding of general relativity 
is incomplete, motivating a number of proposals for modifying gravity at large distances [13-15]. In addition, possible 
alternatives to the dark matter hypothesis have led to theories where gravity changes at slow accelerations or galactic 
scales [16-18]. 



In a previous paper [10] we discussed the possibility of testing general relativity with atom interferometry. We found 
that many relativistic effects will be large enough to be seen in the upcoming generation of experiments. In this work 
we give the details of the framework for calculating the effects of general relativity in an atom interferometer. We 
then apply this to an interferometer in the earth's gravitational field with the motivation of using the high precision 
of atom interferometry to test general relativity in a laboratory experiment. The ability to find GR effects in an atom 
interferometer is more widely applicable as well. In particular we consider other effects such as the Lense-Thirring 
effect. Further, in [19, 20] we apply this technique to find the effect of a gravitational wave. 

We will also discuss a few ideas for strategies to measure several of these GR effects in the lab. We do not attempt 
to prove that such experiments are feasible, since this would require a very detailed analysis of the many relevant 
backgrounds. Instead, we give a few arguments why the most important backgrounds may be controllable. We wish 
mainly to motivate a more careful consideration of these experiments, given the interest in laboratory tests of general 
relativity. 

An attempt was made to make the different sections of this paper as independent as possible. Section II will discuss 
the experimental setup of an atom interferometer and give a description of the usual, nonrelativistic calculation of the 
final phase shift. Section III will describe our method for finding the final phase shift in an atom interferometer in any 
space-time in general relativity. Section IV A will specialize this discussion to the Schwarzschild metric for application 
to an atom interferometer in a weak gravitational field such as the earth's. Section IV B will give our results for the 
GR effects in an atom interferometer near the earth and discuss their physical origin in general relativity. Section 
V will discuss a few ideas for measuring these GR terms in an actual experiment. Section VI will discuss other 
applications of this work including using atom interferometry to measure the Lense-Thirring effect and the expansion 
of the universe. 



Atom interferometry can be used to measure accelerations precisely. In an atom interferometer, the atom is forced 
to follow a superposition of two spatially and temporally separated free-fall paths. When recombined, the resulting 
interference pattern depends on the relative phase accumulated along the two paths. This phase shift is exquisitely 
sensitive to inertial forces present during the interferometer sequence. For example, the goal sensitivity of the next 
generation apparatus exceeds 10~ 14 m/s 2 .[10] 

A single measurement of acceleration in an atom interferometer consists of three steps: atom cloud preparation, 
interferometer pulse sequence, and detection. In the first step, the cold atom cloud is prepared. Using laser cooling and 
evaporative cooling techniques [26], a sub-microkelvin cloud of 10 7 87 Rb atoms is formed in <~ 10 s. Cold atom clouds 
are needed so that as many atoms as possible will travel along the desired trajectory and contribute to the signal. In 
addition, many potential systematic errors (e.g., gravity gradients) are sensitive to the atom's initial conditions, so 
cooling can mitigate these unwanted effects. Here, evaporative cooling is used since it yields the required tight control 
of the initial position and velocity of the cloud that cannot be achieved with laser cooling alone. At the end of the 
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A. Experimental Setup 



3 




FIG. 1: A space-time diagram of a light pulse atom interferometer. The black lines indicate the motion of a single atom. Laser 
light used to manipulate the atom is incident from above (light gray) and below (dark gray) and travels along null geodesies. 
Here the lasers' world-lines are taken to be the two vertical lines on the left and right edges of the graph. 



cooling procedure, the final cloud number density is kept low enough so that atom-atom interactions within the cloud 
are negligible (Section V contains a quantitative discussion of this point). This dilute ensemble of cold atoms is then 
launched with velocity vl by transferring momentum from laser light. To avoid heating the cloud during launch, the 
photon recoil momenta arc transferred to the atoms coherently, and spontaneous emission is avoided [27]. 

In the second phase of the measurement, the atoms follow free-fall trajectories and the interferometry is performed. 
A sequence of laser pulses serve as beamsplitters and mirrors that coherently divide each atom's wavepacket and 
then later recombine it to produce the interference. Figure 1 is a space-time diagram illustrating this process for a 
single atom. The atom beamsplitter is implemented using a stimulated two-photon transition. In this process, laser 
light incident from the right of Fig. 1 with wavevector ki is initially absorbed by the atom. Subsequently, laser 
light with wavevector k2 incident from the left stimulates the emission of a k2-photon from the atom, resulting in 
a net momentum transfer of k c ff = k 2 — ki w 2k 2 . These two-photon atom optics are represented in Fig. 1 by the 
intersection of two counter-propagating photon paths at each interaction node. 

There are several schemes for exchanging momentum between the atoms and the lasers. Figure 2(a) shows the case 
of a Raman transition in which the initial and final states are different internal atomic energy levels. The light fields 
entangle the internal and external degrees of freedom of the atom, resulting in an energy level change and a momentum 
kick. As an alternative to this, it is also possible to use Bragg transitions in which momentum is transferred to the 
atom while the internal atomic energy level stays fixed (see Fig. 3). In both the Raman and Bragg scheme, the two 
lasers are far detuned from the optical transitions, resulting in a negligibly small occupancy of the excited state \e). 
This avoids spontaneous emission from the short-lived excited state. To satisfy the resonance condition for the desired 
two-photon process, the frequency difference between the two lasers is set equal to the atom's recoil kinetic energy 
(Bragg) plus any internal energy shift (Raman). While the laser light is on, the atom undergoes Rabi oscillations 
between states |p) and |p + k e ff) (see Fig. 2(b)). A beam splitter results when the laser pulse time is equal to a 
quarter of a Rabi period (| pulse), and a mirror requires half a Rabi period (n pulse). 

After the initial beamsplitter, the atom is in a superposition of states which differ in velocity by k e fj/m. The 
resulting spatial separation of the halves of the atom is proportional to the interferometer's sensitivity to acceleration 
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FIG. 2: (Color online) Figure 2(a) shows an energy level diagram for a stimulated Raman transition between atomic states |1) 
and 1 2) through a virtual excited state using lasers of frequency u)i and ui2- Figure 2(b) shows the probability that the atom 
is in states |1) and |2) in the presence of these lasers as a function of the time the lasers are on. A -| pulse is a beamsplitter 
since the atom ends up in a superposition of states |1) and |2) while a it pulse is a mirror since the atom's state is changed 
completely. 




FIG. 3: The atomic energy level diagram for a series of sequential two-photon Bragg transitions plotted as energy versus 
momentum. The horizontal lines indicate the states through which the atom is transitioned. The diagonal lines connecting the 
states represent the laser frequencies used in the transition. The result of this transition is to give the atom a large momentum. 



along the direction of k ff. In this work we consider the beamsplitter- mirror-beamsplitter (-| — tt — ^) sequence [28], 
the simplest implementation of an accelerometer and the matter-wave analog of a Mach-Zender interferometer. 

The third and final step of each acceleration measurement is atom detection. At the end of the interferometer 
sequence, each atom is in a superposition of the two output velocity states, as shown by the diverging paths at the 
top of Fig. 1. Since these two states differ in velocity by ~ k e s/m, they spatially separate. After an appropriate drift 
time, the two paths can be separately resolved, and the populations are then measured by absorption imaging. These 
two final velocity states are directly analogous to the two output ports of a Mach-Zender light interferometer after 
the final recombining beamsplitter. As with a light interferometer, the probability that an atom will be found in a 
particular output port depends on the relative phase acquired along the two paths of the atom interferometer. 

Recent atom interferometers have demonstrated sensor noise levels limited only by the quantum projection noise of 
the atoms (atom shot noise) [29]. Assuming a typical time -average atom flux of n = 10 6 atoms/s, the resulting phase 
sensitivity is ~ ^/\/n = 10~ 3 rad/\/Hz and the ultimate phase uncertainly is ~ 10~ 6 rad after 10 6 s of integration. 
This noise performance can potentially be improved by using entangled states instead of uncorrelated atom ensembles 
[30]. For a suitably entangled source, the Heisenberg limit is SNR ~ n, a factor of y/n improvement. For n ~ 10 6 
entangled atoms, the potential sensitivity improvement is 10 3 . Recent progress using these techniques may soon make 
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improvements in signal-to-noise ratio (SNR) on the order of 10 to 100 realistic [31]. Of course, improvements in SNR 
may be easier to achieve simply with increased atom flux without using entanglement. 

Another sensitivity improvement involves the use of more sophisticated atom optics. The acceleration sensitivity of 
the interferometer is proportional to the effective momentum hk e s transferred to the atom during interactions with 
the laser. Both the Bragg and Raman schemes described above rely on a two-photon process for which hk e s — 2hk, 
but large momentum transfer (LMT) beamsplitters with up to lOfifc or perhaps lOOhk are possible [32]. Promising 
LMT beamsplitter candidates include optical lattice manipulations [27], sequences of Raman pulses [33] and adiabatic 
passage methods [34]. Figure 3 illustrates an example of an LMT process consisting of a series of sequential two- 
photon Bragg transitions as may be realized in an optical lattice. As the atom accelerates, the resonance condition is 
maintained by increasing the frequency difference between the lasers. 

Finally, we consider the expected acceleration sensitivity of the coming generation of atom interferometers. We are 
currently building a Rb interferometer that takes advantage of an L ss 10 m vacuum system which allows for a long 
interrogation time of T = 1.34 s [10, 35]. The phase response of the interferometer to an acceleration g is kgT 2 . Using 
10hk LMT beamsplitters and shot-noise limited detection of 10 6 atoms/s for this apparatus results in a sensitivity of 



7 x 10 13 g/>/Hz and a precision of <~ 10 15 g after a day of integration. This estimate is based 



( 10 78Tf^)( 1 -34 s) 2 

on realistic extrapolations from current performance levels, which are at 10 -10 g [36]. 



B. Non-Relativistic Phase Shift Calculation 



In this section we review our non-relativistic method for calculating the phase difference between the two halves 
of the atom at the end of the atom interferometer. These results are well-known [37, 38], but we are not aware of a 
complete, formal derivation of these rules in the literature. Other equivalent formalisms for this calculation do exist 
(see, for example [39, 40]); however, here we review our formalism because the relativistic calculation we discuss in 
Section III has the same structure. For Section III it is necessary to understand the formulae for the phase difference 
(Section II B 1). The proof of these formulae as well as a discussion of their range of validity is given in Section II B 2 
but is not necessary for the rest of the paper. 



1. Phase Shift Formulae 

The main result we will show is that the total phase difference A<p tot between the two paths of an atom interferometer 
may be written as the sum of three easily calculated components: 

A(^tot — ^^propagation ~t~ A(^ S eparation 

+ A(j>i ascr . (1) 

For this calculation and the rest of the paper we take H = c = 1 . 

The propagation phase A(/> pr0 p aga tion arises from the free-fall evolution of the atom between light pulses and is given 

by 

A0 propaga tio„ - ( I {Lc ~ E ^ dt ) ~ S ( / ( ic - E ^ dt ) w 

upper ' lower ^" ll ' 

where the sums are over all the path segments of the upper and lower arms of the interferometer, and L c is the 
classical Lagrangian evaluated along the classical trajectory of each path segment. In addition to the classical action, 
Eq. (2) includes a contribution from the internal atomic energy level E{. The initial and final times t\ and tp for 
each path segment, as well as L c and Ei, all depend on the path segment. 

The laser phase A0i aser comes from the interaction of the atom with the laser field used to manipulate the wave- 
function at each of the beamsplitters and mirrors in the interferometer. At each interaction point, the component of 
the state that changes momentum due to the light acquires the phase of the laser 0z,(io> x c(*o)) = k ' x c(£o) — <^o + <fi 
evaluated at the classical point of the interaction: 

A</w r = ^±0 L (^,x u (^)) - lj2±<f>L(tj,Mtj))) (3) 

V ^ / upper V ^ / lower 

The sums are over all the interaction points at the times tj, and x„(i) and xj(i) are the classical trajectories of the 
upper and lower arm of the interferometer, respectively. The sign of each term depends on whether the atom gains 
(+) or loses (— ) momentum as a result of the interaction. 
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The separation phase A(/> separa ti on arises when the classical trajectories of the two arms of the interferometer do 
not exactly intersect at the final beamsplitter (see Fig. 4). For a separation between the upper and lower arms of 
Ax -X[-x„, the resulting phase shift is 

A</> sep aration = p • Ax (4) 

where p is the average classical canonical momentum of the atom after the final beamsplitter. Of course, the separation 
phase is really an artifact of our semi-classical calculation method based on classical atom trajectories. The observable 
is just the total phase shift at the end of the experiment. 

2. Proof 

The interferometer calculation amounts to solving the Schrodinger equation with the following Hamiltonian: 

Htot =H a + H cxt + Vint (x) (5) 

Here H a is the internal atomic structure Hamiltonian, i? C xt is the Hamiltonian for the atom's external degrees of 
freedom (center of mass position and momentum), and Vint(x) = —fi ■ E(x) is the atom- light interaction, which we 
take to be the electric dipolc Hamiltonian with fi the dipole moment operator. 

The calculation is naturally divided into a series of light pulses during which Vint 7^ 0, and the segments between 
light pulses during which V- ln t = and the atom is in free-fall. When the light is off, the atom's internal and external 
degrees of freedom are decoupled. The internal eigenstates satisfy 

id t \Ai) = H a \Ai) = E t \A t ) (6) 

and we write the solution as \Ai) = \i) e -* B i(*-*o) w jth time-independent eigenstate \i) and energy level E{. 

For the external state solution we initially consider 7J ext = H(x,p) to be an arbitrary function of the external 
position and momentum operators: 

id t |v> = if(x,p)|v>. (7) 

It is now useful to introduce a Galilean transformation operator 

G c = G(x c , Pc , L c ) = e */^ e -ip-x Ce<Pc .i (8) 
which consists of momentum boost by p c , a position translation by x c , and a phase shift. We choose to write 

IV) = G c \4> CM ) • (9) 

We will show that for a large class of relevant Hamiltonians, if x c , p c , and L c are taken to be the classical position, 
momentum and Lagrangian, respectively, then \4>cm) is a wavepacket with (x) = (p) = 0, and the dynamics of \4>cm) 
do not affect the phase shift result (i.e., \4>cm) is the center of mass frame wavefunction) . However, for now we 
maintain generality and just treat x c , p c , and L c as arbitrary functions of time. Combining (7) and (9) results in 

idt\^cM) = {GtH(±,p)G c -iGtd t G c }\cf>cM) (10) 
= {#( x + x c , p + p c ) + p c • x - (p + p c ) • x c + L c } \4> C m) 

where we used the following identities: 

G c xG c = x + x c (11) 

G];pG c = p + p c 

Gjff(x,p)G c = F(x + x c ,p + Pc ) 

Next, we Taylor expand ff (x + x c , p + p c ) about x c and p c , 



H (x + x c , p + p c ) = H(x c , p c ) + VxiJ(x c , p c ) • x + V p i7(x c , p c ) • p + H 2 



(12) 
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where H 2 contains all terms that are second order or higher in x and p. (We will ultimately be allowed to neglect H 2 
in this calculation.) Inserting this expansion and grouping terms yields 



id t 



i>CM) = I (H c - x c • p c + L c ) + (V Xc H c + p c ) ■ x + (V Pc i? c - x c ) • p + i? 2 | \4>cm) 



(13) 



where we have defined the classical Hamiltonian H c = iJ(x c ,p c ). If we now let x c , p c , and L c satisfy Hamilton's 
equations, 



x c = V Pc iJ c 
Pc = -V Xc # c 
L c = x c ■ p c - H c 

with p c = Vi c L c the classical canonical momentum, then \(Pcm) must satisfy 

id t \<i>CM) = H 2 \4>cm) 



(14) 



(15) 



Next we show that it is possible to choose \4>cm) with (x) = (p) = for a certain class of H 2 , so that x c and p c 
completely describe the atom's classical center of mass trajectory. This is known as the semi-classical limit. Starting 
from Ehrcnfcst's theorem for the expectation values of \<Pcm), 



H 2 , Xi 



d^Ho 



d t (xi) = i 



d t { Pi ) = i(\H2,pi }=-{d Xi H 2 



and expanding about (x) and (p), 



d t (xi) 



dt {pi) 



d^Ho 



1 



<x),<P> 



+ dp^^H-2 



(Pj - (Pj)) + d £j d Pi H 2 



+ 



d it H 2 
1 



2! 



<*>,<p> 
dxkdxjdxiH-2 



(x),(P> 

+ d Xi d Xi H 2 



<x>,<P> 

(Pj - (Pj» (Pfc - <Pfc» + 



<x>,<P> 



{Xj - {Xj}) 



(xj - (xj)) + df, 3 d Xt H 2 



(x),(p> 



(x),(p> 

(% - (%)) (xk - (x k )) + 



(x),{p) 



ipj - (pj)) 



(16) 
(17) 

(18) 
(19) 



we find the following: 



d t (xi) = d Pi H 2 
dt (Pi) = - d &i H 2 



+ 7TT dp k dp i dp i H 2 
<x),(p> 2! Pk P] 



<x),<P> 



1 

(x),(p> 2! 



d Xk d Xi d Xi H 2 



<x>,<p) 



Ax% + 



(20) 
(21) 



where Ax^ fe 



(xjXk) — (xj) (£k) and Ap^ fe = (pjpk) — (Pj) {Pk) are measures of the wavepacket's width in phase space 

[72]. This shows that if H 2 contains no terms higher than second order in x and p, then Ehrenfest's theorem reduces 
to Hamilton's equations, and the expectation values follow the classical trajectories. Furthermore, this implies that 
we can choose \4>cm) to be the wavefunction in the atom's rest frame, since (x) = (p) = is a valid solution to Eqs. 
(20) and (21) so long as all derivatives of H 2 higher than second order vanish. In addition, even when this condition is 
not strictly met, it is often possible to ignore the non-classical corrections to the trajectory so long as the phase space 
widths Axjk and Apjk are small compared to the relevant derivatives of H 2 (i.e., the semi-classical approximation). 
For example, such corrections are present for an atom propagating in the non-uniform gravitational field g of the 
earth for which dfdfdfH 2 ~ d^g. Assuming an atom wavepacket width Ax < 1 mm, the deviation from the classical 
trajectory is d t (p) ~ (d^g)Ax 2 < 10~ 20 g, which is a negligibly small correction compared to the <~ 10~ 15 .g effects we 
are considering. 

The complete solution for the external wavefunction requires a solution of Eq. (15) for \(Pcm), but this is non-trivial 
for general H 2 . In the simplified case where H 2 is second order in x and p, the exact expression for the propagator 
is known [41] and may be used to determine the phase acquired by \4>cm) (see also [42]). However, this step is not 
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necessary for our purpose, because for second order external Hamiltonians the operator H2 does not depend on cither 
x c or p c . In this restricted case, the solution for the rest frame wavefunction \4>cm) does not depend on the atom's 
trajectory. Therefore, any additional phase evolution in \4>cm) must be the same for both arms of the interferometer 
and so does not contribute to the phase difference. This argument breaks down for more general H 2 , as does the 
semi-classical description of the atom's motion, but the corrections will depend on the width of \<$>cm) in phase space 
as shown in Eqs. (20) and (21). We ignore all such wavepacket-structure induced phase shifts in this analysis by 
assuming that the relevant moments {Axj k , Apj k} . . .} are sufficiently small so that these corrections can be neglected. 
As shown above for the non-uniform (d 2 g 7^ 0) gravitational field of the earth, this condition is easily met in many 
experimentally relevant situations. 

Finally, we can write the complete solution for the free propagation between the light pulses: 

(x|V, Ai) - (x| G c \cp CM ) \Ai) = j^^e^-^UcMi* ~ x c ) \i) e -«MtF-t,) (22) 

We see that this result takes the form of a traveling wave with de Broglie wavelength set by p c multiplied by an 
envelope function ^cm(x), both of which move along the classical path x c . Also, the wavepacket accumulates a 
propagation phase shift given by the classical action along this path, as well as an additional phase shift arising from 
the internal atomic energy: 

A0 propagation = ( I "( £ c - Ei)dt \ -J2(f F( - L - - E ^ dt ) ( 23 ) 

upper ^Jt' ' lower ^" ' 

where the sums are over all the path segments of the upper and lower arms of the interferometer, and ti, ip, L c , and 
Ei all depend on the path. 

Next, we consider the time evolution while the light is on and Vj n t 7^ 0. In this case, the atom's internal and 
external degrees of freedom are coupled by the electric dipole interaction, so we work in the interaction picture using 
the following state ansatz: 



|*)= /dp^ Ci (p)|V>p>l^) (24) 
J i 

where we have used the momentum space representation of |<^cm) and so \ip p ) = G c e~ 4fl2 ( t_ *°) |p). Inserting this 
state into the Schrodinger equation gives the interaction picture equations, 

id t \V) = ijd P J2 I^p) + H a I*) + ff cxt I*) = H tot |*> (25) 

i 

=* Ci(p) = = \ Jdp'J2 Cj (p') {A t \ (V> P | Mnt(x) IV) l^i) (26) 

where we used (6) and (7) as well as the orthonormality of \Ai) and |V>p)- The interaction matrix clement can be 
further simplified by substituting in \ip p ) and using identity (11): 

^p|Mnt(x)|V P ') - (ple^^VLrt^ + XcJe-^^lp') (27) 
= (p|V-„t(x + x c )|p') e 4 (^^) ( ^ 0) 

where we have made the simplifying approximation that H2 ~ ^ • This approximation works well as long as the light 

pulse time r = t — to is short compared to the time scale associated with the terms dropped from Hi- For example, 
for an atom in the gravitational field of earth, this approximation ignores the contribution m(d r g)x 2 from the gravity 
gradient, which for an atom of size Ax w 1 mm leads to a frequency shift ~ m(d r g)Ax 2 ~ 1 mHz. For a typical pulse 
time t < 1 ms, the resulting errors are < 1 /xrad and can usually be neglected. Generally, in this analysis we will 
assume the short pulse (small r) limit and ignore all effects that depend on the finite length of the light pulse. These 
systematic effects can sometimes be important, but they are calculated elsewhere [43] [44] and do not affect our main 
result for the largest general relativistic effects. 

As mentioned before, we typically use a two photon process for the atom optics (i.e., Raman or Bragg) in order to 
avoid transferring population to the short-lived excited state. However, from the point of view of the current analysis, 
these three- level systems can typically be reduced to effective two- level systems [45] [46] . Since the resulting phase 
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shift rules are identical, we will assume a two-level atom coupled to a single laser frequency to simplify the analysis. 
Assuming a single traveling wave excitation E(x) = E cos (k • x — ut + <p), Eq. (26) becomes 

Ci(p) = \. /rfp'E ^(P'XPl (e^+^-^+V + h.c) |p') e ' (28) 

3 

where the Rabi frequency is defined as fiy = (i\ {—fi ■ E ) \ j) and Wjj = Ei — Ej. Now we insert the identity 

k • (x + x c ) - u>t + (j) = k • x + (k ■ x c (i ) - ut + <pj + £ (k • x c - uj)dt (29) 

boost laser phase Dopplcr shift 

into Eq. (28) and perform the integration over p' using (p| e ±tk x |p') = (p|p' ± k): 

Ci(p) = Skj L(p-k)e i ^e i ^- u+k ^ + ^-^ dt + (30) 

j 

Cj(p + k)e VL e J< ^ J m 2m ' > 

where we define the laser phase at point {io,x c (io)} as (f>L = k • x c (io) — o>to + <P- Finally, we impose the two-level 
constraint (i = 1, 2) and consider the coupling between ci(p) and c 2 (p + k): 

ci(p) = ^c 2 (p + k)e-^ e - j/ * A(p)dt (31) 

c 2 ( P + k) = ^-r»* Cl ( P ) e ^ e j/ * A(p)dt 

Here the detuning is A(p) = loq — lo + k • (x c + JEjL) + the Rabi frequency is O = f2i 2 = (fi 2 i)*, and cj = <*>2i > 0. 
In arriving at Eqs. (31) we made the rotating wave approximation [47], dropping terms that oscillate at (w + w) 
compared to those oscillating at (uj — u). Also, fi^ = since the \Ai) are eigenstates of parity and fi is odd. 
The general solution to (31) is 



ci(p,t) \/ A c (p)e-^ A (P) T -zA s (p) e -^ A (p) T e-^ \ / ci(p,t ) 
c 2 (p + k,i) ) \-iA*(p)ei A (Pl T e^ L A*(p)e3 A (p) T J \ c 2 (p + k, t„) 

A( P ) 



(32) 



A c (p) = cos I - V A(p) 2 + M 2 r I + i sin I - V A(p) 2 + |n|^ r ) (33) 

V2 ^ Va(p) 2 + I^I 2 V2 ^ 

A s (p) = " sinfyA(p) 2 + |»| 2 r) (34) 

^A^ + I^l 2 ^ J 

In integrating (31) we applied the short pulse limit in the sense of k-x c r 2 <C 1, ignoring changes of the atom's velocity 
during the pulse. For an atom falling in the gravitational field of the earth, even for pulse times r ~ 10 fis this term 
is <~ kgr 2 ~ 10~ 2 rad which is non-negligible at our level of required precision. However, for pedagogical reasons we 
ignore this error here. Corrections due to the finite pulse time are calculated elsewhere [48, 49] and they do not affect 
our results to leading order. 

For simplicity, from now on we assume the light pulses are on resonance: A(0) = 0. We also take the short 
pulse limit in the sense of |A(p) — A(0)|r C 1 so that we can ignore all detuning systematics. This condition is 
automatically satisfied experimentally, since only the momentum states that fall within the Doppler width <~ r _1 of 
the pulse will interact efficiently with the light. 



ci(p,t) \_ / A c -iA s e-^\f ci(p,to) A A c = cos 



Mi 

c 2 (p + k,i) J = y -iA* e ^ A c J{ c 2 (p + k,t ) J A s = sin% 
In the case of a beamsplitter (^ pulse), we choose r = f, whereas for a mirror (n pulse) we set r = it: 



(35) 



A , -( 73 vK'M \ ( -ie-* , 
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These matrices encode the rules for the imprinting of the laser's phase on the atom: the component of the atom 
c i(p,io) that gains momentum from the light (absorbs a photon) picks up a phase +</>l, and the component of the 
atom C2(p + k, to) that loses momentum to the light (emits a photon) picks up a phase — fa. Symbolically, 

|p) -> |p + k)e^ (37) 

|p + k) -> |p)e-'* i (38) 

As a result, the total laser phase shift is 

A0i ascr = ^±fefe,x„fe)) - I ^ifofo.x,^)) (39) 

\ ^ / upper \ ^ / lower 

where the sums are over all of the atom-laser interaction points {tj,x u (tj)} and {tj,xi(tj)} along the upper and lower 
arms, respectively, and the sign is determined by Eqs. (37)-(38). 

The final contribution to A0 t ot is the separation phase, A0 sopara ti on . As shown in Fig. 4, this shift arises because 
the endpoints of the two arms of the interferometer need not coincide at the time of the final beamsplitter. To derive 
the expression for separation phase, we write the state of the atom at time t = to + t just after the final beamsplitter 
pulse as 

= |tt u (i)) + |*,(*)> (40) 

where \^ u (t)) and \^i(t)) are the components of the final state that originate from the upper and lower arms, 
respectively. Just before the final beamsplitter pulse is applied, we write the state of each arm as 

|* u (to)) = /dpci(p,to)G u (to)|p)|^i)e w » (41) 

|*,(*o)> = ydpc 2 (p,t )G ; (to)|p)|A 2 )e iei (42) 

where G u = G(x u ,p u , L u ) and G\ = G(x;,p;,L;) are the Galilean transformation operators for the upper and lower 
arm, respectively. These operators translate each wavepacket in phase space to the appropriate position (x„ or x;) 
and momentum (p„ or p;). Here we have assumed for clarity that prior to the final beamsplitter the upper and lower 
arms are in internal states \A\) and |A 2 ) with amplitudes Ci(p, to) and c 2 (p, to), respectively; identical results are 
obtained in the reversed case. We have also explicitly factored out the dynamical phases U and 9i accumulated along 
the upper and lower arms, respectively, which contain by definition all contributions to laser phase and propagation 
phase acquired prior to the final beamsplitter. 

We write the wavefunction components after the beamsplitter in the form of Eq. (24): 

|tt u (i)) = Jd^c[ u \ Vl t)G u \^)\A t ) (43) 

!*«(*)> = /dp£ c fW)G*|P>l^} (44) 

where we invoked the short pulse limit so that e~ lH2T w 1. Next we time evolve the states using Eq. (35) assuming 
a perfect § pulse and using the initial conditions given in Eqs. (41)-(42): namely, c^(p, io) = c i(p> to)e t9u and 
C2^(p,£o) = for the upper arm and c^(p,to) = and Cz\p,to) = c 2 (p, to)e' 9 ' for the lower arm. 

!*«(*)) = |rfp C i(p,to){^G , Jp)|A 1 ) + ^ e ^( x ")G u |p + k)|A 2 )| e ^ (45) 



!*,(*)> = yrfpc 2 (p + k,f )|^e-^ (x!) G ; |p)|^i) + ^G ; |p + k)|yl 2 )| e ^ 



(46) 



We now project into position space and perform the p integrals, 

ci(x-x„,i ) 



(x|* u (t)) = Cl(X Jl u ' oj { e *p»-(*-*») \A,) - ie ^(x.) e i(P„+k)-(x-x„) j e ie u (4?) 
(x|tf,(i)) = C2(X ^|'^ 0) j-z e -^( x <V(P ! - k H x - Xi ) \A X ) + e 4 P''(x-x,) \A 2 )} e ie ' (48) 
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FIG. 4: (Color online) Separation Phase. This is a magnified view of the end of the interferometer which shows the upper and 
lower arms converging at the final beamsplitter at time to, and the resulting interference. The dashed and solid lines designate 
the components of the wavefunction in internal states \Ai) and ^2), respectively. After the beamsplitter, each output port 
consists of a superposition of wavepackets from the upper and lower arm. Any offset Ax = X; — x„ between the centers of the 
wavepacket contributions to a given output port results in a separation phase shift. 



where we identified the Fourier transformed amplitudes using c$(x — x c , to) = Jdp (x — x c |p) Cj(p, to)- The resulting 
interference pattern in position space is therefore 

(x|*(t)) = (x|* u (t)> + <x|*,(t)) (49) 
= ^^{^(x-x^oK^e^-^ 

+ -L \A 2 ) {c 2 (x - x,, t y e <e l P<< x -^ - 1 Cl (x - x u , io ) e ^ e ^(x u ) e4 ( Pll +k).(x-x u) | 

The probability of finding the atom in either output port \A\) or \A%) is 

Kj4i |<x|*(*)>| 2 = |Cl|2 + |C2|2 +^(ic 1 4e' A H C .c.) (50) 

|<A 2 |(x|*(t))| 2 = |Cl|2 + N2 --(icic^ + cc.) (51) 
with ci = ci(x — x„, to) and c 2 = c 2 (x — x;, to). For the total phase shift we find 

A0i = {^ + p„-(x-x tl )}-{^-0 L (x i ) + (Pi-k)-(x-x i )} (52) 

- U - fa - 0z(xi)l + Pi Ax +Ap-(x-x) (53) 



and 



A0p ropaga ti on ,l + A<^i aaori i A0 aopara tion,l 



A0 2 = |fl u + ^l(x„) + (p„ + k) • (x - x„)} - + Pr (x - x ; )} (54) 
, + L (x„)) - 9i + p 2 -Ax +Ap-(x-x) (55) 



A0 prO p a g a tion,2 + A0 laBcri ; 



where pi = p " + ( Pi k - ) and p 2 = ( p " + k)+ Pi are the average momenta in the \A±) (slow) and \A 2 ) (fast) output ports, 
respectively. In general, the separation phase is 

A0 S eparation = P ' Ax (56) 
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which depends on the separation Ax = xj — x„ between the centers of the wavepackets from each arm as well as 
the average canonical momentum p in the output port. We point out that even though the definitions (52) and (54) 
use the same sign convention as our previous expressions for laser (39) and propagation (23) phase in the sense of 
( )u PP or - ( )iower, the separation vector Ax is defined as (x) lowcr - (x) uppor . 

Notice that the phase shift expressions (53) and (55) contain a position dependent piece Ap • (x — x), where 
x = x "+ x ' and Ap = (p„ + k) — p; = p u — (p; — k) , owing to the fact that the contributions from each arm may 
have different momenta after the last beamsplitter. Typically this momentum difference is very small, so the resulting 
phase variation has a wavelength that is large compared to the spatial extent of the wavefunction. Furthermore, this 
effect vanishes completely in the case of spatially averaged detection over a symmetric wavefunction. 

Finally, we show that the total phase shifts A(j>i and A(j>2 for the two output ports are actually equal, as required 
by conservation of probability. According to Eqs. (53) and (55), the contributions to the total phase differ in the 
following ways: 

(A0 prO p aga ti on ,l + A0i ascr ,i) - ^A0 prO p agat ion,2 + A</>i asCTi2 ) = <j>L (xj ) - (j) L (x„ ) (57) 

= k- (xj - x„) = k- Ax 



A(/> separ ation,l - A(/> scparat i on ,2 = Pi ' Ax - p 2 ■ Ax = -k • Ax (58) 

Together these results imply that A^i = A0 2 and prove that the total interferometer phase shift A0 tot is independent 
of the output port. 

The accuracy of the above formalism is dependent on the applicability of the aforementioned stationary phase 
approximation as well as the short pulse limit. The stationary phase approximation breaks down when the external 
Hamiltonian varies rapidly compared to the phase space width of the atom wavepacket. The short pulse limit requires 
that the atom's velocity not change appreciably during the duration of the atom-light interaction. Both approximations 
are justified to a large degree for a typical light pulse atom interferometer, but in the most extreme high precision 
applications such as we consider here, important corrections are present. However, we emphasize that these errors due 
to finite pulse duration and wavepacket size are well-known, previously established backgrounds. Although they must 
be accounted for experimentally, these corrections do not affect the leading order general relativistic effects which we 
seek to calculate. 



III. GENERAL RELATIVISTIC DESCRIPTION OF ATOMIC INTERFEROMETRY 

We are interested in the leading order effects of general relativity in an atom interferometer. Section 11 B described 
our non-relativistic calculation. Here we build on the results of that Section to create a formalism which treats the 
entire calculation in a relativistic manner. It is very difficult to solve for the quantum mechanical evolution of the 
atom in a general metric background. Thus, we will use the semi-classical approximation method outlined for the 
non-relativistic calculation. This method can be used in relativity, with some minor modifications, since in general 
relativity the concept of potential is replaced by the least action principle. 

In brief, the method is as follows. Using the prescription in Section II B 1, the free propagation of the atoms and the 
light in Figure 1 is treated non-quantum mechanically. Thus, both the laser pulses and the atoms are taken to move 
along geodesies of the space-time. The description of the atom-light interaction is taken from non-relativistic quantum 
mechanics, but must be described in a covariant manner as will be discussed below. Finally, the total resulting phase 
shift must be a coordinate invariant. Further, to understand the result it is necessary to write it in terms of the 
physical variables of the problem as measured by an experiment, thus removing any coordinate dependence from the 
answer. 

Our objective is to calculate the leading order general relativistic effects in order to to explore interesting and 
possibly measurable signals. In order to simplify the calculation, many sub-leading order effects will be dropped 
including effects due to the finite pulse time of the lasers, AC Stark shifts, and the errors in the semi-classical 
approximation due to the finite size of the atom's wavefunction. All these may give important backgrounds but they 
can be and have been calculated easily in the non-relativistic formalism. We are interested in the largest effects that 
arise due to general relativity and so we can ignore all these effects. 

We will consider an atom interferometer in a background space-time with metric g M „. The proper time for a particle 
in this space-time is then given by dr 2 = g^dx^dx" . We will take h = c = 1. 
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A. Dynamics of the Interferometer 



The trajectories of the atoms and the laser pulses are found by solving the geodesic equation 

d 2 x^ _ dx a dx< 3 „ 

^ + r ^7" ' (59) 

where T is the affine connection and Greek indices run to 3. In order to compute the leading order GR effects we will 
calculate the phase shift using the approximations explained in Section II, for which it is sufficient to find the motion 
of the center of the atomic wavefunction. Of course, there are corrections to this semi-classical (or stationary phase) 
approximation due to the finite size of the atomic wavefunction which have been discussed in Section II B. However, 
the leading order GR effects are just large enough to be experimentally measured and therefore these corrections 
to our GR results are negligible and will be ignored. The interferometer is then defined by the initial space-time 
position and velocity of the atom before the first beamsplitter pulse, and by the starting positions of the three laser 
pulses which define the interferometer, the beamsplittcr-mirror-beamsplitter sequence. Once these are given, the rest 
of the interferometer is found by calculating the intersection of the geodesies as shown in Fig. 1. We assume atomic 
transitions only occur when the atom is simultaneously in the presence of both laser beams. We will assume that for 
every atom-light interaction point, the right laser is always turned on sufficiently far before the left laser so that the 
atom is already in the presence of the light from the right laser when the light from the left laser hits it. Thus the 
atomic transitions will always occur when the light from the left laser reaches the atom. This choice makes a small 
but potentially measurable effect that will be discussed below. We will call the left laser the 'control laser', because 
its timing determines the timing of the interferometer and the right laser the 'passive laser'. 

The intersection of the initial beamsplitter pulse (from the control laser) and the initial atomic trajectory defines 
point A. After the first beamsplitter interaction, the half of the atom which is not affected by the light travels on 
to B along the same trajectory. The half of the atom which is affected travels on a new trajectory originating from 
point A but with a new velocity which is a function of the incoming momenta of the atom and the light pulse as will 
be described below. The two halves then travel to B and C respectively, which are defined by the intersections with 
the mirror pulse. At B, the atom gets a kick from the light in the same way as before. At C, it loses momentum 
via stimulated emission. These halves then travel on to D and E respectively, where they interact with the final 
beamsplitter and interference is assumed to occur as described above. The relativistic calculation of the final phase 
shift is described below. Thus, given the initial conditions for the atom and laser pulses, the interferometer can be 
calculated in a fully covariant manner. 

We must now give a coordinate invariant description of these initial conditions. This means they must be written 
in terms of the physical variables measured by the experimentalist. For all the calculations described here we will 
consider the laser to be at a fixed spatial coordinate location, xl- As we will see, for all the metrics we consider, we 
will choose coordinates such that this is a suitable model for the laser. The initial beamsplitter pulse then defines the 
start of the interferometer and should be considered to be at an arbitrary time coordinate, t\ (unless the experimenter 
somehow has independent knowledge of the metric). The mirror pulse is then taken to leave the laser at xl a time 
T later as measured on the laser's clock. Similarly, a time T later on the laser's clock the final beamsplitter pulse is 
emitted from the laser. The laser's proper time is given by 

T= fdr= f y/g^dt = f y/g^dt. (60) 

Solving these two equations for t 2 and t 3 yields the time-coordinates at which the mirror and final beamsplitter pulses 
originate from the laser. This then defines the three laser pulses in a coordinate invariant way The atom's initial 
position may be changed depending on the application and thus it is harder to give a single, relativistic description of 
it. One natural way to define it is to take the atom to begin at the same position as the laser (and thus at a well-defined 
time on the laser's clock), and then to travel for a certain amount of time (again on the laser's clock) before the first 
beamsplitter pulse is emitted. This defines the initial position in a relativistically invariant manner. The atom's 
initial velocity can also be defined in several ways. For example, it could be taken as the radar ranging velocity that 
the laser sees. Where it is relevant, we will usually consider it to be the more experimentally realistic velocity that 
would result from getting some large number of momentum kicks from the laser. Once the initial conditions for the 
atom interferometer have been defined in a coordinate invariant manner, we can calculate the entire interferometer 
sequence. 

The atom- light interaction can most easily be thought about in a local Lorentz frame, x' , (LLF) of the atom 
(essentially Riemann normal coordinates). This is a choice of coordinates such that the atom is at rest at the origin 
of these coordinates and space-time is locally flat near the atom. Specifically in the LLF the metric is locally flat 
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with vanishing first derivatives at the position of the atom so that 

= (61) 



1 a' /3' 



x'=0 



and near the origin (the position of the atom) 



g^','(x')= r)„'w+0{x' 2 ), (62) 

where r\ is the flat metric. In these coordinates, the leading order effects of the interaction with the light on the atom 
are just the non-relativistic quantum mechanical rules given in Section II. So in this frame the spatial momentum of 
the atom after the transition is equal to the sum of the spatial momenta of the atom before the transition and the 
incoming light which causes the transition, namely 



At' 1 



dr 



m fi 



dT 1 

ux atom 



dT 
after 



ax light 



d\ 

before 



(in LLF) (63) 



where 2 = 1,2, 3, m m and mg are the rest masses of the atom before and after the atomic transition, and A is the affine 
parameter for the laser's null geodesic (the analogue of r in Eq. (59) but with different units). Typically m; n = nifi±uj a 
where w a is the frequency difference between the initial and final (internal) atomic states that are coupled by the laser 
interaction [73]. Note that in the case of a multi-photon transition (e.g., the two-photon transitions shown in Fig. 1) 

the photon momentum ^ ht in Eq. (63) must be replaced by the sum of the momenta of all photons that contribute 
to the transition. Since the metric deviates slightly from flat, there are small corrections due to the tidal effects of 
gravity over the size Ax of the atomic wavefunction which are 0(RAx 2 ), where R represents the Riemann curvature 
tensor. Thus any effects of gravity are suppressed by the size of the atomic wavefunction, so even the leading order 
Newtonian contributions to Riemann can be neglected during the atom-light interaction (see the discussion of the 
semi-classical approximation in Sec. II B 2). For a description of transforming coordinates to the LLF see [50]. 

To find these spatial momenta in the LLF (x 1 ) it is necessary to know them in the main coordinate system (x) 
being used and then transform them to the LLF. The atom's spatial momenta are simply determined from its geodesic 
and proper time coordinate r. However, a light ray travels on the same geodesic no matter what its momentum is. 

We must then know the correct affine parameterization, a^ ght (A), of the null geodesic such that x ^'| ht is actually the 
spatial momentum and not just proportional to it. This is determined by the initial conditions for the laser pulse. 

However, we cannot simply take the spatial momentum dx ^ ht of the laser pulse at emission to be equal to k l , 

XL 

the laser's wavevector, because the coordinates x l do not necessarily have any physical meaning and so such a choice 
would be coordinate dependent. We must instead write our answers only in terms of physical variables; in this case 
we must use the k which an observer would measure the laser to have. For definiteness we will assume the laser's 
frequency is measured by an observer at the same position xl as the laser and not moving with respect to it. Then 
this defines 



= cu = k (64) 



dx" 

where u> is the frequency of the laser and — d ° ba is the four-velocity of the observer. This equation implicitly 
defines the A such that our observer sitting on the laser (in the laser's LLF) measures the emitted pulses to have the 
proper frequency uu. This pulse is then propagated to the atom in the main coordinate system of the problem (x). 
When the light reaches the atom, its momentum is transformed to the atom's LLF (x 1 ) in order to find the momentum 
transferred to the atom. Eq. (64) is critical as it modifies the result for the GR corrections at leading order as we 
will see explicitly below. 

Note that this rule (Eq. (63)) of adding the momenta only applies in the LLF, and there is in fact some ambiguity 
about which LLF to use: the rest frame of the atom before or after the interaction. This ambiguity implies relativistic 
corrections to the recoil velocity that are 0(v 2 ), thus changing the atom's velocity at O(vf). This is far too small to 
be visible in the experiments we are considering; however if it becomes necessary to compute such corrections, the 
formalism given here would have to be adapted. 

We have now completely determined the dynamics of the interferometer in a fully relativistic framework. Im- 
portantly, the initial conditions have been given in terms of physical variables of the problem and so there is no 
coordinate dependence left. All that remains is to determine the rules for calculating the final phase difference in 
general relativity. 
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B. Relativistic Phase Shift Formulae 



The propagation phase is proportional to the integral of the Lagrangian, i.e., the action. In general relativity, the 
action for a particle moving in a background space-time is the length of its world-line, S = J mdr, where the mass m 
of the particle is irrelevant for the equations of motion but is inserted so that the Lagrangian reduces to the normal 
non-relativistic Lagrangian in the appropriate limit. Since S = J Ldt, the Lagrangian is L = mj^. To demonstrate 
that this reduces to the expected Lagrangian in a weak gravitational field, insert the Schwarzschild metric (Eq. (71)) 

and take a post-Newtonian expansion in <fr and v 1 = ^g-. Then L = m — j-mv 2 + iruf> + 0(v 4 , cf) 2 , v 2 <p) as expected. 
By analogy with the non-relativistic formula in Section II, the general relativistic action is the propagation phase 



^propagation 



j Ldt = j mdr = j p^dx 1 * (65) 



where p^ = m^- is the particle's momentum. The last equality follows from p^dx* 1 = rng liv ^ r dx v — m^dr. Notice 
that this is the opposite sign convention from the non-relativistic expression for propagation phase from Section II. 

The separation phase follows essentially from the formula outlined in Section II applied in the LLF. The separation 
phase in the LLF is EAt' — p ■ Ax 1 . We then employ the standard trick, to write the formula in this frame as a 
coordinate invariant 



^separation 



D 



E 



p^dx" 1 ~EAt'-p-Ax* (66) 



where the integral is taken along the null geodesic connecting points E and D (the final beamsplitter pulse from laser 
1). Since this is a coordinate invariant and is true in the Local Lorentz Frame (x'), it is then valid in all frames. Then 
in any coordinate system (x) the separation phase is 



(•D 

^separation 



IE 

Here p M is the average momentum of the two halves of the atom at points D and E 



/ %dx?- (67) 

J E 



-u 1 ( dx^ 

P =2 ™°^ 



dx» 
+ m — — 
dr 



(68) 



and the momenta are evaluated at points D and E after the final beamsplitter pulse and in the same output port 
(either slow, \A\), or fast, \A%), whichever is being used to compute the final phase shift; see Fig. 4). Here m is the 
mass of the atom in the chosen output port. The small coordinate ambiguities in this formula are negligible for every 
metric considered, as will be shown below. As is clear from the formulas, the separation phase (67) can be thought 
of as the last piece of propagation phase (65). The term EAt is then roughly the phase acquired by bringing the 
half of the atom that transitioned earlier at the second beamsplitter up to the same time as the second half of the 
atom, so they can interfere. We have chosen to define the separation phase along the null geodesic which is the final 
beamsplitter pulse because this is the first hypersurface on which interference can be considered to have occurred. Of 
course, since quantum mechanics is linear, we can choose any later hypersurface, add the two halves of the atom's 
wavefunction and calculate a separation phase then (of course we would also need the correct rule for propagating 
the atoms forward after the final beamsplitter). In other words, in order to find the total phase shift, the atom's 
wavefunction can be evaluated at any point in space-time after the final beamsplitter pulse. For example, the number 
of atoms in each output port is often measured with a detection laser pulse. This would correspond to a final null 
geodesic which is equivalently good for calculating the final phase shift. 

As shown in Section II, the laser phase imparted to the atom during a beamsplitter or mirror pulse is the phase 
of the light at the interaction point. Since this applies in the LLF and is also a coordinate invariant statement, it 
can be applied in any frame. A null geodesic is a line of constant phase for the laser since e.g. it is the world-line 
in space-time that a crest of the laser pulse follows. Thus the light's phase at an interaction point is the same as its 
phase at the time of emission of that light pulse from the laser. For example, the laser phase of the pulse from laser 
1, the control laser, at point A is just the phase of laser 1 at time t\ in Fig. 1. The total laser phase shift is then 
the sum of the laser phases from each laser over all the interaction points. For example, in the slow output port, the 
total laser phase is 

A</Or = Mh) - 20ii(t 2 ) + - 0L 2 (4 2) ) + 0L2(4 2) ) + fe(4 2) ) - 0L2(4 2) ) (69) 
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where (j>Li(t) means the phase of laser i at time t. Because we have assumed that laser 2 is always on [74], the 
transition points are defined by the intersection of the atom's geodesic with the light pulse from laser 1. Of course, 
the interaction takes some small amount of time, but we assume that the phase of the laser at the beginning of the 
interaction is the one imparted because the corrections to this approximation are typically small and do not affect 
the leading order GR result (see the discussion of the short pulse limit in Section II B 2). Thus the phase which laser 
2 imprints on the atom is the phase of the light from laser 2 (the passive laser) which is passing the atom at the first 
instant of the interaction as set by laser 1 (the control laser). Finding this phase requires determining the time that 

(2) (2) 

this part of the light left laser 2, which we label t A , t B , etc. We will usually assume that the phase of laser i is just 
its frequency (as defined by Eq. (64)) times its proper time <pLi{t) = ^iTLi(t)- For a time independent metric this 
is also 4>Li{t) = ^iy/goot. So in this case, the contributions to laser phase from laser 1, the control laser, completely 
cancel. 

We now have rules for finding the final phase shift in an atom interferometer, 

A</>tot = A0 prO p aga tion + A0 scpara ti on + A0i ascr (70) 

in general relativity. 

It can be seen that this formula is independent of the output port used to calculate the phase. The propagation 
phase does not depend on output port. Laser phase depends on output port since in the slow output port it includes 
the phase of the laser at point D, while in the fast output port it includes the phase of the laser at point E. As can be 

seen from Figure 1, this is a difference of <~ ui (t^ — t^^j where ui is the frequency of laser 2. The separation phase 

also differs between ports because the momentum used in Eq. (67) is the average momentum of the relevant component 
of the atom's state after the final beamsplitter pulse. In the fast output port this momentum is ~ mv r « k e fi larger 
than in the slow output port. So the difference in separation phase between the fast and the slow output port is 
~ k e ff(xE — %d) ~ 2lo(xe ~ %d), since k e s « 2lo for a two-photon transition. But xe — xe ~ Ie — since points 
D and E lie on a null geodesic, and as a result we find that xe — xe ~ \ (t^ — t^^j . Therefore, the difference in 

the separation phase between the two output ports is exactly compensated by the change in the laser phase, and Eq. 
(70) is independent of which output port is used. 

IV. GR EFFECTS IN THE EARTH'S GRAVITATIONAL FIELD 

The methods of the previous section can be used to find the effects of general relativity in an atom interferometer 
in the earth's gravitational field. The space-time can be modeled with the Schwarzschild metric 

ds 2 = (1 + 20) dt 2 - -\-dr 2 - r 2 dn 2 (71) 

where — is the gravitational potential. For simplicity, in this section the rotation of the earth is neglected. 

It will not modify the GR effects given here, and the possibility of measuring relativistic effects associated with 
that rotation will be considered in Section VI B. Of course, this rotation will contribute non-relativistic terms that 
can be backgrounds, which can easily be calculated using simpler, non-relativistic methods. Because the earth's 
gravitational field is weak, ~ 10 -9 at the surface, we can take a post-Newtonian expansion in 0. In order to study 
effects beyond GR as well, we will calculate the phase shift in the parameterized post-Newtonian (PPN) expansion 
of the Schwarzschild metric 

ds 2 = (1 + 20 + 2f3<j) 2 )dt 2 - (1 - 2 7 0)dr 2 - r 2 dfl 2 . (72) 

Here (3 and 7 parameterize modifications of general relativity and j3 = 7 = 1 gives normal GR. For the results 
presented in this section we will generally be considering this PPN expansion of the Schwarzschild metric, though 
there is little difference in a weak gravitational field. Section IV A contains the details of the phase shift calculation 
for this metric. In Section IV B we present the results of this calculation and explain their physical significance. It 
can be read without the previous Section IV A. 

A. Interferometer Calculation in the Schwarzschild Metric 

The geodesic equation (59) for the metric (72) cannot in general be solved exactly. We approximate the solution 
using a power series solution in t. By varying the order of the series, we ensure that we use a sufficiently high order to 
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include all measurably large terms in the final phase shift. For simplicity, we present first the approximate solutions 
using the metric (71) for a radial geodesic: 

r(r) = r + v r0 r - | (0 r 0(r„)) r 2 - ^v r0 (5 2 </>(r )) r 3 + G(r 4 ) (73) 

, . . V^ro + ^ + ^Po) Vro^v^o + V + 2r]^(r Q )d r ^(r ) 2 
t(r) = t H : , „ J?- \ t 7 1 , KT7Z ^ + (74) 



l + 20(r o ) ' (l + 20(r o )) 

^ o + ?? + 2 ^ (ro) ((4 v 2 o + rj + 2 ^(r )) d r ^{r f - « r 2 (l + 2^(r o ))a 2 0(r o )) 



r 3 



3(l + 20(r o )) 3 



+ G(r 4 ) 



where r\ = <7/^^f-^f- is for null geodesies and 1 for time-like geodesies. Note that these are roughly the normal 
parabolic paths with some relativistic corrections. Also, light rays do 'bend' under gravity, but in these coordinates 
that effect shows up in the equation for t(r) only. In these solutions, the potential <j) has effectively been expanded 
around ro, making this approximation better the closer the paths are to ro or, roughly, the smaller r is. The scale 
this is to be compared to is the radius of the earth, Re, as this determines the size of the higher r derivatives of <f>. 
Since the atom travels a distance much smaller than Re, this approximation works very well here. The paths in the 
PPN metric (72) are 

, (-^ 2 o (-l+7)+^ + 2(^ 2 o (/3-2 7 ) + /3r ? )0(r o )-6t> r V7^ro) 2 )^(r o ) 2 

r { r) r + v r0 r+ 2(-l + 27^ro))(l + 20(r o ) + 2/30(ro) 2 ) 

1 ( 27 K 2 o(-l + 27) - V - 2 K 2 q(/3 - 3 7 ) + 0rj) + 8v 2 r0 f3^(r ) 2 ) d r </>(r ) 2 
6 Vr0 [ (1 - 2 7 ^(r )) 2 (1 + 2^(r ) + 2/3<^(r ) 2 ) 

2 (-Vro -V + 2^ 2 7^M) (-2 + /3 + 7 - 6(/3 - 7)^0) - W ' 37)0M 2 + 16/3 2 7 ^(r ) 3 ) c^M 2 | 

(l-27^r o ))2(l + 20(r o ) + 2/30(r o )2) 2 ~ ~~~ + 

v 2 rQl dMro) (1 + -WM) (% 2 o + g ~ 2^ 2 o 70Po)) 5 2 0(ro) \ 3,-,,^ ^ 
1 - 2 7 ^(r ) + (-1 + 2 7 #r„)) (1 + 2<^(r ) + 2/3</>(r ) 2 ) ) ^ 

. , , jjg^Q "ro(l + 2^(ro))y ^- a g^ff a ^(ro) 2 

"" ' " ' W l + 20(ro) + 2/3^(r o ) 2T 1 + 20M + 2/3</>M 2 T 




2 ^° 70(ro) ((-^ 2 o(-5 + 2/3 + 7) + V + 2 K 2 (-67 + 0(8 + 7)) + Wv) 0M + 



l + 2^(r o ) + 2/30(ro) 2 
2/3 ( W 2 (8/3 - 19 7 ) + 2/3?7) 0(r o ) 2 - 36f 2 o /3 2 7 0(r o ) 3 ) d r </>(r ) 2 + 
v 2 r0 (-l + 2 7 0(r o )) (1 + 2(1 + /?)«Kr ) + 6/3^(r ) 2 + 4/3 2 </>(r ) 3 ) #0(r„))) r 3 )/ 

(3 ((-1 + 2 7 0(r o )) (1 + 2#r„) + 2/30(r o ) 2 ) 2 ) ) + 0(r 4 ) (76) 

These geodesies are calculated on a computer with a symbolic algebra package and so can easily be found to higher 
orders. We present here the results up to third order to illustrate the method without overcomplicating the equations. 

Theoretically, the intersections points of the laser pulses with the atom geodesies can now be found as discussed 
above. However, these geodesies are too complicated to solve exactly so we must make two approximations. First, the 
equations are solved self-consistently only to the order in r to which the entire calculation is done. In practice this 
means, for example, inverting the series to find r(i) which can then be plugged in to x(t) to give x(t) for the atom. 
An analogous procedure is used to find x(t) for the light. The atom and light coordinate trajectories are then equated 
and solved perturbatively order by order in t. Second, we must still expand in the variables which are small in order 
to simplify the results sufficiently so that they remain tractable. The correct way to do this expansion that keeps only 
relevant terms and avoids an "order of limits" problem is to Taylor expand in all the variables simultaneously, taking 
into account their relative sizes. Specifically, let e signify 10 _1 and plug in e 9 </> for <fi everywhere (since <j> <~ 10~ 9 ), 
and similarly for the other dimensionless variables: v L ~ 1(T 7 , %f ~ v r ~ lCT 10 , ^ ~ lCT 15 , ^ ~ 1(T 2 , and 
mT ~ 10 26 . Since the corrections to the parabolic paths are small, we have a few large terms at low orders which 
give easily solved equations and then very many small corrections. So we can Taylor expand all our results in e and 
keep only terms which are large enough to possibly affect the final answer. For the intersection points, we usually 
keep any terms which are > 10 -30 i?£. We vary this order to make sure we haven't neglected any relevant terms in 
the final phase shift. The entire calculation is done on a computer using a symbolic algebra package. 
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We take each laser to be at a fixed coordinate location because the lasers are assumed to be fixed to the earth. 
Note that this is not a geodesic. A fixed radial coordinate implies a fixed physical position since this is a static, 
time-independent metric. Effects such as time variations in the earth's gravitational field or oscillations of the laser 
platform which cause the laser's effective coordinate location (or the whole metric) to vary are very small and so 
will affect the leading order GR signal at an unmeasurable level. Of course, depending on the phenomenological 
characteristics of the signal being searched for, such effects crossed into the zeroth order Newtonian signal could be 
relevant backgrounds. They can then be calculated simply in a non-relativistic fashion. For now we ignore them since 
we are interested in calculating the effects of general relativity. In Section V we consider measuring these GR signals 
and there we discuss the relevant backgrounds. 

The rest of the calculation will be illustrated using metric (72). The interferometer is defined, as shown above, by 
the initial positions and momenta of the laser pulses and the atom. For simplicity, and because it illustrates all the 
effects we will be interested in, we make the following choices. The lasers will be at positions r = and r^ 2 with 
frequencies u>\ and u>z, with laser 2 above laser 1. We will take the usual definitions k e fi = oji +1^2 and u c g = u>\ — LU2- 
The pulses from laser 1 will originate at ti = 0, ti = , ===== and t 3 = 2t 2 . Following the prescription 

y/l+2<f>(rL 1 )+2(34>(r Ll ) 2 

given in Eq. (64) for referencing the light momenta, an observer at the laser has coordinate velocity 

jjp, = dx obs = ( ^obs dx ohs \ = ( 1 \ 

dr \ dr ' dr J [ ^TTW^J+WW^f ' / 
giving an initial light momentum for the pulses from laser 1 of 

^ff - (^f > ) = «i (l + <Krm) + l<t>(r Ll f mr Ll f, 1 - 7^) + fy 2 *^) 2 ) . (78) 

This is close to (u>i,u>i) but with small GR corrections. These corrections will modify the GR effects in the final 
answer at leading order and so they must be included. 

We will take the atom to be initially at r = tl ± at t = 0. For now we leave all expressions in terms of the unphysical 
coordinate launch velocity, vl = fjr- We show below that this makes no important difference to the final phase shift 
and it keeps the expressions simpler. The geodesies and intersections can now all be found as explained in Section 
III. 

As one example, we describe the calculation of intersection point C. The initial velocity of the fast half of the atom 
at point A is found by adding the light momentum to the atom's initial momentum as described above. For the 
atom-light interaction at A, the coordinate transformations to the LLF take a vector to 

V& F = W'V (79) 

with 

b ( > = y/(l + vl + 2vl^{r Ll )) (1 - 2^{r Ll ) + 2(3^(r Ll f) (80) 

bl = -v L {l + 2^{r Ll )) (81) 

bl = vl V(1 + gggj) (1 ~ 20(r £l ) + 2ggj ~^j (82) 

b\ = + 27^)) (l + v 2 L + 2vh<l>(r Ll )) (83) 

The velocity of the half of the atom going from A to C in the main coordinate system (metric (72)) is then 



dr 
dr~ 



k e f[ k cii v L k cS 2 k cS uj a uj cS fc off 3fc cff 2 2 

v L + + -^o- + v L v L -yd r <p(r Ll ) + 7 d r <p(r Ll ) (84) 

to 2m z to to to mm 2 to 



which is roughly vl + v r with relativistic corrections. Intersecting the atom and light geodesies using the approxima- 
tions described above gives point C. The expressions for the coordinates are very long and so are given in Appendix 
A. Of course, these are roughly just rc ~ + (vl + ) T and tc ~ T as they would be in the non-relativistic 
case. 

(2) (2) 

Computing the laser phase requires finding the times t A , t B , etc (see (69)) by dropping null geodesies from these 
intersection points to the second laser at Tl 2 - 

The entire calculation is done on a computer using a symbolic algebra package (Mathematica) so all Taylor series 
orders, initial conditions and such can be changed easily. The calculation was written for a general metric so the same 
code is used to calculate the phase shift for several different metrics including (71) and (72). 
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B. General Relativistic Effects and Interpretation 

We present here the results of the calculation of the phase shift in an atom interferometer placed in a weak 
gravitational field such as the earth's. The physical origins of the important terms in the phase shift will be discussed, 
focusing on the GR terms and their interpretation. 

To summarize the variables we are using, T is the interrogation time between pulses on the laser's clock, Wj is the 
frequency of laser i, k tf = u>i + w 2 , w e fj = U\ — lu 2 , co a is the frequency of the atomic transition between states 1 and 
2 of the atom (see Figure 2(a)), m is the rest mass of the atom in the lower ground state, Vl is the atom's launch 
velocity in the radial (vertical) direction, is the position of laser i, <j>(r) is the local gravitational potential (for a 
perfect Schwarzschild metric (f>(r) = —^-), g = — V</>, and 7 are PPN parameters in the metric, and c = h = 1. 
We present the results for the final phase shift for metric (72) in Table I. The phase shift has been expanded into a 
sum of terms and we have grouped terms that have the same scalings with experimental control parameters. Table 
I also displays the results of a non-relativistic (NR) calculation for the phase shift in a gravitational potential (see 
Section II for a description of the NR calculation) . Note that we have kept the laser frequencies the same for all three 
pulses, though in a real experiment these would have to be tuned to keep the transitions on resonance in order to 
compensate for the Doppler shift due to the atom's acceleration under gravity [75]. 

To understand the GR effects underlying some of these phases, recall that, roughly, the atom interferometer is 
sensitive to accelerations. The following discussion will be highly coordinate dependent and not rigorous, but its only 
purpose is to gain some intuition for the GR effects we find. Combining the geodesic equations (59) for the spatial x 1 
(i = 1, 2, 3) and t, the coordinate acceleration of an atom in the frame of Eq. (72) is 

dv — * — ' — * 

— = -VW + (8 + 7)0 2 ) + 7(3(tf • ff - 2v 2 )V<p + 2v{v • V</>) (85) 
at 

with v = ^| for this equation only. The acceleration is approximately that from Newtonian gravity, —V<j>, but with 
leading order GR corrections. These corrections fall into two classes, both of which will interest us. 

The V^ 2 terms are related to the non-linear nature of gravity, the fact that a gravitational field seems to source 
itself in general relativity. This could also be called the non-Abelian nature of gravity since gravitational energy 
gravitates through a three-graviton vertex. To see that this is the origin of the V0 2 terms, note that, because of 
these terms, the divergence of the gravitational field given in Eq. (85) is nonzero (here by gravitational field we mean 
g = 4| from Eq. (85)). Just as for an electric field, a nonzero divergence of the gravitational field implies a local 
source density (in general relativity this means a local energy density) that is proportional to that divergence. So Eq. 
(85) implies that there is a local energy density in free space proportional to V • g oc V • V0 2 = 2(V</>) 2 . But note 
that to leading order V</> ss g so that V • g oc g 2 . In other words, the local energy density is proportional to the field 
squared, exactly as expected from the electric field analogy. This energy is then the source of the V0 2 terms. The 
non-linearity of gravity is parametrized in the standard way by the PPN parameter j3. 

The other terms in Eq. (85) proportional to iPVcf) are velocity dependent forces. These terms are related to the 
gravitation of the atom's kinetic energy (or the kinetic energy of the source mass in the frame where the atom is 
stationary and the source is moving), since all energy, not just mass, gravitates in general relativity. 

The non-linear GR corrections in Eq. (85) are smaller than Newtonian gravity by a factor of </> <~ 1CT 9 , while the 
velocity dependent force terms are smaller by v 2 <~ 10~ 15 for the atom velocities we are considering. We will see 
that the non-linear terms can only be measured through a gradient of the force produced and so are reduced by an 
additional factor of fi 10m ss 10~ 6 for a 10m long experiment. Both effects are then ~ 10 _15 g. 

These effects can be seen in the total phase shift in the interferometer. Table I presents the answer for the total 
phase shift as found by the relativistic calculation outlined above. It lists all the terms in the total phase shift large 
enough to be measured by the initial apparatus. Effectively, the local gravitational acceleration is expressed as a 
Taylor series in the height above the Earth's surface. The first phase shift in Table I represents the effect of the 
leading order (constant) piece of the local acceleration while the 2nd and 10th terms are the next gradients in the 
Taylor series. Notice that even the second gradient of the gravitational field is relevant for this interferometer. The 
terms in this list that have been measured agree with the results of previous experiments. The largest two phase 
shifts due to the first two terms in the Taylor expansion of the local g field were known and measured several times 
(e.g. [8, 28]). The 4th and 5th terms arise from the second order Doppler shift of the laser's frequency as seen by the 
moving atom. These Doppler shift and finite speed of light corrections (terms 4 and 5) were known and measured 
to cancel each other in a 'symmetric' interferometer [51]. The 7th term is proportional to the two-photon detuning 
between the difference in the two lasers' frequencies, w ff, and the resonant frequency of the atomic transition, ui a - In 
any practical experiment this detuning is kept quite small and this term will be negligible. The terms proportional 
to u> a and Lo c g almost cancel since the lasers' frequencies are usually chosen to be on resonance with the atomic 
transition and so could never have been measured given the precision of previous interferometers. These terms were 
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GR Phase Shift 


Size (rad) 


Interpretation 


NR Phase Shift 


1. 


-k cS gT 2 


3. x 10 s 


Newtonian gravity 


-k cS gT 2 


2. 


-k e s{d r g)v L T 3 


-2. x 10 3 


1st gradient 


-k cB {d r g)v L T 3 


3. 


-^k eS (d r g)gT 4 


9. x 10 2 




-^k eS (d r g)gT 4 


4. 


-3k ctt g 2 T 3 


-4. x 10 1 


finite speed of light and 




5. 


~3k ctt gv L T 2 


4. x 10 1 


Doppler shift corrections 




6. 


2 


-7. x 10 _1 


1st gradient recoil 


2 

-iS(9r 9 )T 3 


7. 


(UJ CB - UJa)gT 2 


-4. x 10" 1 


detuning 




8. 


(2 - 2/3 - y)k cH g<j>T 2 


-2. x 10" 1 


GR (non-linearity) 




9. 


- 3 4?9T 2 


2. x 10" 2 






10. 


-^k aB v 2 (d 2 g)T 4 


8. x 10 -3 


2nd gradient 


-^k cB v 2 (d 2 g)T 4 


11. 


-fk cfi (d r g)gv L T 4 


6. x 10" 4 






12. 


-4k eS (d r g)v 2 L T 3 


-3. x 10" 4 






13. 


2co a g 2 T 3 


2. x 10 -4 






14. 


2tu a gv L T 2 


-2. x 10 -4 






15. 




7. x 10" 6 


2nd gradient recoil 


- 7 -^v L (d 2 g)T 4 


16. 


-12k cB g 2 v L T 3 


-7. x 10" 6 






17. 


-7k cS g 3 T 4 


4. x 10 -6 






18. 


-5k cH gv 2 L T 2 


3. x 10 -6 


GR (velocity-dependent force) 




19. 


(2-2/3- 7 )fc cff 9 r (#)w L T 3 


2. x 10" 6 


GR 1st gradient 




20. 


^(4-4/3-3 7 )fccff0(9rff)0T 4 


-2. x 10 -6 


GR 




21. 


(w e ff — wj (d r g)vLT A 


2. x 10 -6 






22. 


(ueff - w ) {d r g)gT 4 


-1. x 10" 6 






23. 


-^(2-2/3- 7 )fc cff <? 3 T 4 


-3. x 10" 7 


GR 




24. 


-^ {dr g)v L T 3 


-2. x 10" 7 






25. 


- 2 -^(d r g)gT 4 


2. x 10 -7 






26. 


m " 


-1. x 10" 7 






27. 


6(2 - 2/3 - >y)k c s4>g 2 T 3 


5. x 10 -8 


GR 




28. 


3(uj cB -uj a )g 2 T 3 


4. x 10 -8 






29. 


3 (Weff — Wo) gVhT 2 


-4. x 10" 8 






30. 


6(1 - f3)k cB( f>gv L T 2 


3. x 10" 8 


GR 





TABLE I: A list of all the terms above a certain size in the phase shift from the full GR calculation for metric (72), along 
with their numerical size in radians and an interpretation. The NR phase shift column shows the results of a completely non- 
relativistic phase shift calculation for comparison. The sizes of the terms assume the initial design, sensitive to accelerations 
~ 10~ 15 s, which has L = 9 m, T = 1.3 s, v L = 13 f , k cH = 2j^^, co a = 6.8 GHz, and m = 81 GeV (for 87 Rb). All detuning 
terms assume lu c h — u) a — 1 kHz. Note that there is some ambiguity in how some of the terms are grouped since by definition 
g — —d r <f>- 



not previously known because their calculation requires a fully relativistic calculation. The recoil shift -^T d r g was 
known and measured [52, 53]. 

The 8th, 18th through 20th, 23rd, 27th, and 30th terms arise only from GR and are not present in the results of our 
Newtonian calculation. The 8th and 19th terms arise in part from the non-linear nature of gravity. This is clear since 
they look like the analogue of kgT 2 and the 1st gradient terms but with g replaced by the part of the acceleration 
coming from the non-linearity of gravity, g<j) in Eq. (85). Similarly the 18th term arises in part from the velocity 
dependent forces in Eq. (85) since it appears to be an acceleration ~ gv\. Note that of course the acceleration from 
these velocity dependent forces is actually proportional to the integrated effects along the entire trajectory of the 
atom. However, this is obscured by the expansion we are taking, and so we just see the largest term, oc v\, with terms 
proportional to the other velocities v r , gT and so on, farther down the list. In fact, for every term oc Vl, we expect 
and see a term with vl replaced by gT also on the list, since the velocity changes over the course of the interferometer 
by roughly this amount. 

We now address the fact that we have left the phase shift in Table I in terms of the unphysical (coordinate- 
dependent) launch velocity. This is the only coordinate dependent variable in Table I; all others (fc e ff, T, etc.) have 
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Parameter 


Total Phase Shift 


Propagation Phase 


Separation Phase 


Laser Phase 






Dependence 


Coefficient 


Coefficient 


Coefficient 


Coefficient 


Size (rad) 


1. 


k ctt T 3 {d r g) 





1 


-1 





4 x 10 10 


2. 


kcttgT 2 


-1 


-1 


1 


-1 


3 x 10 8 


3. 




1 


2 


-2 


1 


3 x 10 3 


4. 


CO a gT 2 


-1 


-1 








3 x 10 3 


5. 


k eS (d r g)T 3 v L 


-1 


2 


-2 


-1 


2 x 10 3 


6. 


k ce {d r g)<t>T z 





27 + 2/3-2 


-2 7 -2/3 + 2 





3 x 10 1 


7. 


k eB gT 2 v L 


-3 


-5 


5 


-3 


1 x 10 1 


8. 


k e sg4>T 2 


2-2/3-7 


2-2/3-7 


-2 + 2/3 + 7 


2-2/3-7 


2 x 1CT 1 


9. 


k cS g 2 T 3 v L 


-12 


-17 


17 


-12 


7 x 1(T 6 


10. 


k cB dr{gct>)T 3 v L 


2-2/3-7 


-4 + 4/3 + 2 7 


4-4/3-27 


2-2/3-7 


2 x 1(T 6 


11. 


k cB gT 2 vl 


-5 


-9 


9 


-5 


5 x 10 -7 



TABLE II: A breakdown of some of the terms in the phase shift, Table I, list by origin. The sizes are given for the coefficient 
of the term in the total phase shift. Note that there is some ambiguity in dividing terms 9 and 10 since by definition g = — <9 r </>. 



coordinate invariant definitions. Ultimately, in any real experiment the experimenter determines how to measure the 
launch velocity, and this gives the physical, coordinate-invariant definition. On this point, different experiments will 
surely vary, so here we assume a simple prescription but leave Table I in the general form in terms of Vl, which 
should allow any prescription to be applied. We assume that the atom is launched by the lasers using til photon kicks 
(Raman or Bragg transitions). The experimental definition of the physical launch velocity will then be v p = njj^-. 
Note that til, k e s, and m all have physical, coordinate- indepedent definitions. Repeating the normal atom- light 
interaction calculation (Eq. (63) and ensuing discussion) til times then gives the relation 

vl « v p ( 1 - ^ - ^ (1 - 7 - r ±-l^L (i _ 1(t>) \ (86 ) 
\ m 2 m J 

with higher order terms dropped. It is not surprising that there are higher order GR corrections when the coordinate 
launch velocity is written in terms of a physically measurable parameter. Here we are only interested in this if it 
changes the GR effects we seek to measure, for example by changing the dependence on the PPN parameters. It 
is clear that substituting Eq. (86) into the phase shift in Table I will not affect the two GR terms we are most 
interested in, 8 and 18. It can affect other GR terms at the level of 10~ 6 rad and below, but it cannot remove 
totally the dependence on the PPN parameters of GR terms 8 and 18, and so it does not qualitatively change their 
interpretation. 

The origin of some of the largest terms in the phase shift list highlights important differences between a relativistic 
calculation and a non-relativistic one. From Table II we see that frequently the contributions to a given term in 
the phase shift lift from the propagation and separation phases cancel, and so the term can be considered to come 
from laser phase. There is even a term in the propagation and separation phases that is larger than kgT 2 which 
cancels. Note that this term is k c gcT 3 (d r g) if we do not take c = 1. Many of the terms in propagation and separation 
phase can be considered to arise from the fact that, in a relativistic calculation, the endpoints of the interferometer 
(points D and E in Fig. 1) are not simultaneous. This has a large effect because the phase of the atom evolves at 
a rate proportional to its mass m, as follows from the separation phase formula (67). In other words, the Compton 
wavelength measures separation in time just as the de Broglic wavelength does in space. 

It is interesting to consider an atom interferometer with only a single laser driving the atomic transitions directly 
between two levels, instead of the normal two-photon transition through a virtual intermediate level. In this case, 
there is no passive laser and the laser phase is zero, as mentioned in Section 111 B. This means that the k e ggT 2 term 
is removed as are most of the k e g terms. However we must now have the two atomic levels spaced by an energy which 
is roughly the frequency of a laser, so ui a ~ k, instead of u e s. After doing such a single laser calculation we find that 
the largest term in the phase shift is u a gT 2 , which is roughly the same size as kgT 2 . This phase shift now comes 
from propagation phase instead of laser phase, as would be guessed from Table II. This term arises because the rest 
mass of the atom is different in the two atomic states. This means that the dominant phase shift does not depend on 
the laser frequency and is instead set by the intrinsic structure of the atom. In the normal non-relativistic calculation 
there is only the k e sgT 2 term which comes from laser phase (not from propagation phase), and there is no co a gT 2 
term at all. So the non-relativistic single laser calculation gets the origin of the major part of the phase shift wrong, 
and it gives an answer which is off by the amount the laser frequency u> is detuned from resonance uj a . 
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The term -5k cS gT 2 vl (term 18) receives contributions from the velocity-dependent forces in Eq. (85), but its 
coefficient is independent of 7. There are two canceling contributions to this term coming from the 7 terms in the 
force law for the atom and the photon. Note that the 7 terms in the equation of motion for the light are not suppressed 
by any small velocity factors (since c = 1), so they are just as large as the 'normal' Newtonian gravity term (this 
is the origin of the famous factor of 2 in the equation for the bending of light by the sun). This phase shift term 
thus measures both the effect of gravity on light and the velocity-dependent force on the atom. If we put a different 
parameter, <5n g ht, in front of the </> in the component goo of the metric governing the motion of the light and redo the 
entire calculation, this term becomes (4 + <5n g ht + 7n g ht — 1atom)k c sgT 2 vf j7 where the 7's are the PPN parameters in 
the metrics for the light and the atom. This term then tests a matter-light principle of equivalence, namely that they 
both feel the same metric. A new force may well couple to light and matter in a different manner than gravity and 
so may be testable in this way. 

V. MEASUREMENT STRATEGIES 

In this section we consider several ideas for measuring GR or beyond GR effects using atom intcrferomctry. Possible 
measurement strategies are discussed for testing the Principle of Equivalence (PoE) and measuring the effects of the 
velocity dependent forces arising from GR, the non-linear nature of gravity, and the bending of light in a gravitational 
field. We have discussed some of these ideas in a previous paper [10] but here we give a few more thoughts. Our 
main motivation is to describe some phenomenological differences between the GR effects we have calculated and 
classical, non-relativistic effects in order to show that these GR effects are not coordinate artifacts, and are possibly 
distinguishable from Newtonian gravity and other backgrounds. We do not claim to have proven that every conceivable 
background is under control. Instead we only wish to argue that it may be possible to test GR in the lab using this 
technology, and so it is well worth trying to design experiments to do so. 

In general we consider a~10m long interferometer with T ~ 1 s and <~ 10 6 atoms cooled and launched per shot, 
resulting in a final phase sensitivity of ~ 10~ 6 rad after about 10 6 shots. This implies sensitivity above atom shot noise 
of 10~ 15 g where g is the acceleration due to gravity on earth. As we have seen, this is enough to start measuring GR 
effects in the lab. Of course, there are many possible improvements to this technology, and the potentially achievable 
sensitivities (at least above shot noise) were discussed in [10]. 

A. Principle of Equivalence 

Atomic interferometers have been used to test the Principle of Equivalence (PoE) [8, 58]. We intend to perform 
a similar test to 10 -15 accuracy using a Rb apparatus currently under construction (see Section II) [35]. The basic 
idea is to co- locate two atomic clouds of different species or isotope (e.g., 85 Rb and 87 Rb) and run simultaneous 
atom interferometers using the same laser pulses on these two clouds (a simultaneous dual species fountain has been 
demonstrated [59]). This test is essentially looking for the presence of new forces that act in a PoE violating way, 
namely that are not proportional simply to the mass of the atoms but instead depend on the atomic species or isotope. 
It seems likely that any new force would be PoE violating since it is hard to design a force other than gravity that 
does not violate the PoE. Indeed most examples of new forces in the literature do violate the PoE. This experiment, 
being done on the surface of the earth, has maximal sensitivity to PoE violating forces with a range greater than the 
earth's radius. An infinite range force would also be detectable in this manner. As the range, A, of the new force 
decreases below the earth's radius, the sensitivity of the experiment decreases as a ~ 10" 15 (^) (with the overly 
simplistic assumption of uniform earth density) since the volume sourcing the new Yukawa force goes as <~ A 3 . Notice 
that this does give sensitivity to forces with shorter ranges, down to ~ 10 m, the rough size of the experiment and 
the distance over which it is not clear what the local mass distribution is. Figure 5 shows projected limits from this 
experiment on a new force that couples to Baryon number (B) and Baryon minus Lcpton number (B - L), and arising 
from a light dilaton [62]. The large suppression factors come from the fact that 85 Rb and 87 Rb have similar couplings 
to these new forces and also roughly equal masses. Also shown in Figure 5 are the current experimental limits on 
a force coupling to B coming from an equivalence principle experiment [61]. This line shows the difference that a 
realistic earth model makes for the limits. 

One possible background to the PoE measurement arises from the fact that the earth's gravitational field is not 
constant but in fact has a gradient ~ where R e is the radius of the earth. Thus an increase in height of 5h ~ 10 nm 

changes the acceleration by <~ g^- ~ 10~ 15 g, which is at the desired sensitivity. The primary way to reduce this is 
to use a differential measurement between two simultaneous atom interferometers run with the same laser pulses. In 
the PoE measurement it is necessary to co-locate the two isotopes in the initial atomic trap to < 10 nm in order to 
avoid this systematic. The differential sag in the combined gravity plus magnetic trap depends only on the isotope 
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FIG. 5: (Color online) Projected limits to be set by the Principle of Equivalence measurement. Limits are shown (labeled solid 
lines) for several possible new Yukawa forces arising from a coupling to Baryon number (B), Baryon minus Lepton number 
(B-L) and the dilaton. These assume a uniform earth density. Previous limits are shown in solid (yellow) shading [60]. The 
dashed line is the current limit on a force coupling to B from an Equivalence Principle experiment using a realistic earth model 
[61]. 



masses since we can trap Rb and Rb using states with the same magnetic moment. For our designed magnetic 
trap curvature, the expected differential sag is ~ 10 nm. Further, this differential position can be measured by 
characterizing the magnetic trap in situ, which can be used to reduce the systematic phase error below 10~ g. 



B. Velocity Dependent Forces 



In order to measure the GR effects of the 'velocity dependent forces' and the gravitational effect on the laser light, 
we must pick out the phase shift term with the v\ scaling. The basic idea is that it is very difficult for classical gravity 
to mimic the effect of a true velocity dependent force. Again we will consider a differential measurement between 
simultaneous atom interferometers since this aids greatly in the control of many backgrounds. Here we consider two 
atom interferometers with differing launch velocities in order to select the velocity dependent term we are looking 
for. This will naturally cancel many phase shifts, including the leading order one from gravity, k c sgT 2 . Of course, 
there are still several possible background terms that are larger than the signal and scale with launch velocity. These 
include terms coming from the earth's gravity gradient, k g(d r g)T 3 VL, the effect of the doppler shift of the laser as 
seen by the moving atom, k c ggT 2 VL, and the effect of the Earth's rotation, k c gVL^ 2 T 3 , where is the component of 
the Earth's rotation rate perpendicular to the launch velocity [38] . Magnetic fields do not induce velocity dependent 
phase shifts when the internal atomic state remains unchanged, as in a Bragg beamsplitter. All of these background 
terms, although much larger than our signal, scale differently with vl and T than our signal. There are no terms 
which scale as v 2 L T 2 \ this is a unique sign of GR. Varying these around the typical values (vl = 13— and T = 1.3 s) 
then allows the GR term to be picked out from the backgrounds with a sensitivity limited only by the atom shot noise. 
It is crucial for this fit that Vl can be known experimentally very precisely (better than the ratio of the background 
to the signal). This is possible since the launch velocity is precisely linked to laser frequencies (see the discussion of 
physical velocity and Eq. (86)) which can be known extremely well. 
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Additionally, in case backgrounds do become a problem, it is possible to reduce the measured size of the background 
terms even before this fit. Because the GR term scales as (w^T) 2 , there is no loss in the signal by going to the regime 
where the launch velocity is large. In this regime, the atom's velocity is roughly constant over the length L of the 
interferometer and L — vlT. If vl and T are then always scaled inversely so that L is fixed, then the signal does not 
change but all the background terms do. Further, by taking vl large, all the background terms are suppressed by at 
least one power of vl because they all have more powers of T than of vl- For example, the gravity gradient term 
becomes k e s(d r g)^. Thus, the sizes of the phase shifts from these backgrounds can be directly reduced even before 
data analysis. 

C. Non-Linearity of Gravity 

To measure the terms which arise from the non-linearity of gravity, it may help to run three simultaneous atomic 
gradiometers along three mutually orthogonal axes in a 'divergence configuration'. Such a configuration effectively 
measures the divergence of the local gravitational field, which must be zero in Newtonian gravity outside the source 
mass. This should then allow the non-linear GR effect, k^g^T 2 , to be picked out. In particular, the atoms can 
be launched along a single, large (e.g. 10 m) vertical axis. Then three perpendicular atom interferometers can 
be run along this same axis using three perpendicular sets of lasers. Thus the atoms can be split vertically or in 
either horizontal direction to make the three perpendicular atom interferometers. Yet all three interferometers are in 
essentially the same position, separated only by the much smaller v r T ~ 1 cm. 

One question which arises is the extent to which the three laser axes can be made mutually orthogonal, since they 
must be perpendicular to one part in 10 9 (since on the earth <f> ~ 10~ 9 ) in order to reduce the earth's gravity gradient 
below the GR signal. Methods for measuring angles with nanoradian precision have already been demonstrated, albeit 
for angles much smaller than 90 degrees [54] . It may also be possible to geometrically construct laser axes which are 
perpendicular to very high accuracy by using the ability to accurately measure distances with a laser interferometer. 

It is also possible to ameliorate the requirement on the perpendicularity of the laser axes by reducing the local gravity 
gradient with an appropriately constructed local mass distribution. We will show that it is possible to reduce the 
gravity gradient along all three perpendicular axes by 0(1) of its natural size on the earth. With such a construction, 
the three gravity gradients can then be measured by the atom interferometer itself. It should then be possible to 
make minor modifications to the mass distribution to cancel the gravity gradients with increasing precision, without 
an exact knowledge of the angles of the atom interferometers. Every order of magnitude reduction in the size of 
the local gravity gradient reduces the requirement on perpendicularity of the laser axes by an order of magnitude. 
Since the atom interferometer can measure gravity gradients very precisely, it may be possible to align the lasers to 
sufficient accuracy without a complicated alignment mechanism. 

Now we must show that the earth's gravity gradients in all three directions (i.e. d x g x , d y g y , and d z g z ) can be 
cancelled to 0(1). If the z-axis runs perpendicularly to the local surface of the earth at the point in question, then 
g z ~ — G M carth [761 is negative and the gradient d z g z w —2 „ g * is positive. Adding more mass on the z-axis either 

"earth -"earth 

completely above or completely below the atom interferometer will only add a positive quantity to d z g z , thus increasing 
the gravity gradient. Therefore we must add mass around the atom interferometer. As a proof of principle, take a 
sphere of mass centered on the point in question with a cylindrical hole along the z-axis (the atom interferometer 
apparatus would be placed in this hole). Assuming the radius of the cylinder is small compared with the radius of 
the sphere, the Newtonian gravitational field due to the sphere inside itself has a derivative with an opposite sign 
compared to the earth's gradient: d z g z = — |7rGp sp hcrc- Since in vacuum the divergence is zero in Newtonian gravity 
{d x g x + d v g v +d z g z = 0) and there is a rotational symmetry about the z-axis, the other two components of the sphere's 
gravity gradient are d x g x — d y g y — —\d z g z . With an appropriate choice of /9 S phcro then, the earth's gravity gradient 
can be cancelled off. Notice that the sphere cancelled the z-component of the earth's gradient while the cylindrical 
hole cancelled the x- and y-components, since in the infinite limit the cylinder only has x- and y-components of 
acceleration. 



D. General Backgrounds 

One possible background to any measurement made with an atom interferometer arises from the interaction of the 
atoms with ambient electromagnetic fields. The electric fields present in any realistic setup are too small to give 
detectable phase shifts as they are easily screened. Surface effects such as the Casimir interaction are negligible as the 
atoms are kept far from all surfaces. Only ambient magnetic fields can give large enough phase shifts to be potential 
backgrounds. An atom responds strongly to a background magnetic field, so we usually consider a magnetically 
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FIG. 6: (Color online) The response function, T gz , as defined in Eq. (87). It gives the phase shift response of the atom 
interferometer to a Fourier component of the local gravitational field with wavelength A and amplitude 10 -9 gr oa rth- All curves 
assume the example 10 m atom interferometer. The top curve assumes it is run in a 'symmetric' configuration (see text), the 
next lower curve assumes the atoms are dropped from rest at the top of the device and the following two curves assume a 
downward launch velocity of 1— and 10— respectively. 



shielded interferometer. Ambient magnetic fields can be shielded down to the nT level (see for example [55]) which 
leads to appropriately small phase shifts. The atoms are prepared in a magnetically insensitive, m = 0, state and so 
the energy shifts that arise are second order in the magnetic field, AE = \ctB 2 where a is the second order Zeeman 
coefficient of the state. Since the internal levels of the atom can have different values of a, magnetic phase shifts 
are generally smaller when the internal state is fixed, as is the case in an interferometer that uses Bragg transitions. 
Systematic phase shifts can then only result from spatial variations in the field: 

A(f) sa ~k cS — B Q — — T 
m oz 

where B — Bq + SB, Bo ~ 100 nT is the constant bias magnetic field and SB is a small field perturbation. Variations 
of SB ~ 1 nT over the length of the interferometer give negligible phase shifts for a 87 Rb apparatus. The above 
formula generally holds for field perturbations that vary on length scales that are long compared to the interferometer 
arm splitting. In the opposite limit, the interferometer averages over perturbations with wavelengths that are small 
compared to the overall interferometer length. This spatial averaging behavior also occurs for short wavelength gravity 
perturbations, as we describe in more detail below. 

Another potential background in the interferometer arises from atom-atom collisions within the atom cloud. For 
87 Rb, the frequency shift of the atomic state \F = 2, mp = 0) due to atom-atom collisions is 



Sv k. (-0.9 mHz) 



10 9 cm" 3 / V 1 MK 



for a cloud of number density n and temperature r [56, 57]. Unlike in atomic clocks, phase errors due to this effect in 
a f _ 7T _ § pulse sequence atom interferometer implemented using Bragg atom optics are the result of unequal densities 
between the two arms. Nominally, the upper and lower arms of the interferometer have the same atom number density, 
but an imperfect initial beamsplitter can cause an asymmetry between the arms, resulting in a phase shift 
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for a density difference An. This represents an upper bound on the atom-atom phase shift, since in reality the cloud 
density decreases in time during the experiment due to ballistic expansion [77]. Making the conservative assumption 
that the density difference can be controlled at the level of ^ <~ 1 CP 2 implies a phase error of ~ f 0~ 4 rad. However, 
this systematic offset is not a concern for many of the experiments we consider since it is expected to cancel as a result 
of our differential measurement strategies. This cancellation relies on the condition that the density does not depend 
on any of the other control parameters in the experiment, an assumption that must be verified experimentally [78]. 

An additional tool that may be useful for reducing backgrounds and picking out the signal terms is the possibility 
of running the atom interferometer with different pulse sequences. While the §-7r-§ sequence is an accelerometer, 
more complicated sequences can be used which suppress accelerations and leave the gravity gradients, or vice versa. 
In general, either T 2 or T 3 terms, or both, can be removed by a suitable choice of pulse sequence [40, 63]. This can 
remove most of the relevant backgrounds to the velocity-dependent force measurement since they scale with higher 
powers of T. These different pulse sequences improve the ability to pick out a term that scales in a particular way 
with the control parameters. 

Finally, it is possible that small gravitational 'anomalies' due to local masses may be a background. In fact, only 
the small wavelength variations in the local g field can be a relevant background for these experiments. Perturbations 
to local g at wavelengths larger than the rough size of the experiment (e.g. 10 m) are well described by the Taylor 
series expansion that we have assumed for the earth's field. Since these perturbations are naturally small compared 
to the earth's field, the differential measurement strategies discussed above for the Principle of Equivalence, velocity 
dependent force, and non-linearity of gravity measurements will remove these long wavelength perturbations in exactly 
the same way as they removed the earth's field. 

Short wavelength perturbations in the local gravitational field can be relevant. Luckily, the atom interferometer 
naturally averages over these perturbations because of its spatial length. We write the total phase shift due to gravity 
anomalies along the vertical (z) direction summed over all wavelengths A as 



where the gravitational field at a position z in the interferometer is Sg z (z) — J 5g z (X)e^'dX, Sg z (X) is the Fourier 
component of a gravity perturbation with wavelength A, and T gz {X) is the interferometer's gravity perturbation 
response function. Figure 6 shows the response of an atom interferometer in the example 10 m configuration to 
perturbations in the local g field of wavelength A. The top curve is for the atom interferometer run in the 'symmetric' 
configuration where the atoms are launched from the bottom of the interferometer region with exactly the right 
velocity to stop at the top, Vl = gT and the first and last beamsplitter pulses happen when the atom is at the bottom 
on its way up and on its way down, respectively. The next lower curve assumes the atom is dropped from rest at 
the top. The next two curves assume the atom is launched downwards from the top with velocities 1™ and 10^, 
respectively. As expected, the atom interferometer always averages down the perturbations on scales below its size, 
here 10m. If the atom is launched so that its velocity is never zero then the suppression is much bigger at shorter 
lengths. For the lower two curves, the phase falls off as A 2 , as compared with A for the case in which the atom comes 
to rest during the interferometer. Roughly, the more uniform the atoms' velocity, the larger the reduction that comes 
from averaging over the small scale gravity perturbations. This is a very useful tool for reducing backgrounds from 
local masses. It also scales favorably with the length of the interferometer, since a longer interferometer suppresses 
larger distance scales, leading to an even greater suppression at short distances. 

The differential measurement strategies suppress the longer wavelength contributions. In fact, at very long wave- 
lengths there is no difference between what we have called a 'perturbation' and the previously included part of earth's 
gravitational field. The differential measurement strategies were designed to allow us to control systematics arising 
from the earth's gravitational field. For the proposed Principle of Equivalence measurement, atom clouds of the 
two isotopes are separated by less than l^iin. This provides a large suppression to longer wavelength gravitational 
perturbations. Similarly, in the divergence strategy measurement (see Section V C) the three atom interferometers are 
separated by ~ 1 cm. For wavelengths longer than this scale the divergence measurement looks like a true divergence 
(instead of a finite difference approximation to the derivative) and so the longer wavelength gravitational modes are 
suppressed. As discussed above, the velocity dependent force measurement benefits greatly from the ability to launch 
the atoms with a large initial velocity to suppress all Newtonian gravitational effects without suppressing the GR 
signal. This suppression of all wavelengths can be seen by comparing the lowest two curves in Figure 6. The lowest 
curve has a larger launch velocity and is therefore suppressed at all wavelengths as expected. 

If we know the magnitude of the perturbation in the local gravitational field, we can compute the phase shift induced 
in an interferometer from Eq. (87). This implies some constraints on the mass and distance of nearby objects that 
must be taken into account when designing an actual experiment. Fig. 7 displays Sg z (X) for several example sources. 
Each source is taken to be some perpendicular distance from the center of the atom interferometer. We have always 
plotted the magnitude of the Fourier transforms in Figures 6 and 7 (i.e. we have taken the sum in quadrature of the 
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FIG. 7: (Color online) The magnitude of the power spectra of the local gravitational field, Sg z {\) from Eq. (87), for several 
example sources. The solid (blue) curve is a 10~ 2 kg point source, 10 cm from the center of the interferometer. Similarly, the 
dotted (purple) curve is a 1 kg source at 1 m and the dash-dotted (yellow) curve is 1000 kg at 10 m. The long-dashed (green) 
curve is a thin 10 m long rod of mass 10 kg, parallel to the interferometer, whose center is 1 m from the interferometer. 



sine and cosine components). An object that is farther away than its size looks like a point mass. Its gravitational 
field has a Fourier transform, 8g z (X), that is sharply peaked around its distance from the interferometer. An object 
that is large compared to its distance (represented by the rod in Figure 7) produces a Sg z (X) that is sharply peaked 
around the object's size. This analysis allows us to predict the expected size of local gravitational perturbations as a 
function of wavelength. 

The differential response curve allows us to compute systematic errors arising from the specific gravity environment 
of our interferometer. Quantitative estimates of these effects requires knowledge of the local dg z (z), which may be 
obtained through a combination of modelling and characterization. The atom interferometer itself can be used as a 
precision gravimeter for mapping 5g z {z) in situ. By varying the interrogation time T, a local gravity measurement 
can be made over a small spatial region. Many such measurements in different positions can be made by varying the 
launch velocity and time of the initial beamsplitter, resulting in a measurement of 8g z (z). 

The interferometer's response to the short wavelength modes can be quite suppressed as in Fig. 6. Additionally, 
the amplitudes of the gravitational perturbations can be kept small as in Fig. 7. Further, depending, to some extent, 
on the nature of the gravitational source, the phase shifts due to these local gravitational anomalies arc unlikely to 
vary at order one with the control parameters used in the experiment, or to vary at order one from shot to shot. 
This leads to a further suppression, since a truly constant phase shift would not be a background to many of these 
proposed experiments. Thus, it seems possible to reduce the background phase shifts due to local gravitational effects 
below the required 10 -6 rad level. 

There are, of course, many details we have not addressed here that are important for a real experiment. Here we 
have given our ideas for ways of distinguishing the main GR effects of the earth's gravitational field from the relevant 
backgrounds. These effects are found to be possibly phenomenologically distinguishable from Newtonian gravity and 
other backgrounds, and are thus true GR (or beyond GR) effects, indescribable in a Newtonian model. We have tried 
to motivate why these or similar experiments may be possible, but we have not proven that they will work in the full 
detail necessary for a real experiment. We just wish to motivate future work. 



28 



VI. OTHER GR EFFECTS? 
A. Hubble Expansion (Vacuum Energy) 

The cosmological expansion of the universe can be described by the metric 

ds 2 = dt 2 - e 2Ht dx 2 (88) 
where H is the Hubble constant. Using the geodesic equation (59), a particle moving in this metric has an acceleration 

d 2 x dx dt , s 

For laboratory atomic velocities this would be <~ 10 _18 g, large enough to be very interesting compared to the initial 
sensitivity of the interferometer under construction ~ 10~ 15 g. We would like to consider whether it is then possible 
to measure the expansion rate of the universe in an atom interferometer. Since the local neighborhood has collapsed 
gravitationally, the normal cosmological expansion of the metric is presumably not occurring inside the solar system. 
However, a true cosmological constant would necessarily be present everywhere and must therefore affect the metric 
locally. We will loosely model the local effects of the vacuum energy using metric (88) (to include the effects of local 
mass a McVittie metric [65] should really be used), but it must be remembered that H will then presumably refer 
only to the contribution of the vacuum energy. 

Unfortunately, this is just an artificial, coordinate acceleration, not the physical acceleration that would be measured 
by any experiment. A calculation of the acceleration that would be measured in an actual experiment, either an 
observer using radar-ranging to determine the atom's position or using an atom interferometer, shows that the 
physically observable accelerations are of 0{H 2 ). 

Another way to see an effect linear in H is to consider two observers on geodesies of (88) which are at constant x 
positions. The spatial distance (defined by integrating only the spatial terms in the metric at a fixed time) between 
these observers is e Ht Ax and so it would appear that there is a relative velocity of O(HAx). In fact a radar-ranging 
calculation of the distance or velocity shows essentially the same thing. This is not a surprise, since this is in fact just 
the normal observation of Hubble's law for the expansion of the universe. Such observers represent galaxes, which are 
measured to recede from each other at velocities of order H. However we cannot create a similar laboratory version 
of this experiment in order to measure the Hubble constant (really vacuum energy) without starting the atom with 
an initial velocity of O(H), which defeats the purpose. The galaxies naturally have such velocities, and so tend to 
mark the cosmological expansion, due to the action of the 0(H 2 ) acceleration acting for the age of the universe H -1 . 

More generally, it seems there cannot be effects linear in H, at least in 'normal' variables such as position, velocity, or 
acceleration, in any experiment we can set up. Without access to observers in some special reference frame such as the 
galaxies provide, any experiment will be free-falling in the cosmological expansion. By the Principle of Equivalence, 
such an experiment will then see a metric that is locally flat with corrections proportional to the Reimann curvature. 
But this is R cx H 2 (all components are either 0, — H 2 , or ±H 2 e 2Ht ). This is very similar to transforming Eq. (88) 
to the static patch 

ds 2 = (1 - H 2 r 2 ) dt 2 1 dr 2 -r 2 dn 2 . (90) 

In these coordinates it is clear that all accelerations (all Christoffel symbols), the Reimann curvature and other similar 
quantities will be 0(H 2 ). Of course, the argument that in a LLF the answer will just be proportional to the curvature 
and therefore to H 2 misses the important point that the Reimann tensor (at least in either of the coordinates (88) 
or (90)) still has coordinate dependence, and therefore cannot actually be observable. This problem highlights the 
importance of working entirely in physical variables, as we stressed for the atom interferometer in Section III. 

The expansion of the universe is occasionally proposed as an explanation for the Pioneer Anomaly [64] which is an 
anomalous acceleration of the Pioneer spacecrafts of order H. Radar ranging the Pioneer is quite similar in spirit to 
the atom interferometer, which can be thought of as laser ranging the atom to find its acceleration. Thus, similar 
arguments apply to this case as well, and the Pioneer anomaly cannot be explained within general relativity as being 
due to the cosmological expansion. 

Intuitively, it is impossible to observe effects linear in Hubble in a local experiment because of the equivalence 
principle. Essentially, everything in the experiment is 'falling' together in the expansion of the universe. Similarly, it 
is impossible to detect (at leading order) the acceleration toward the dark matter of the galaxy. Only the gradient of 
this force is detectable, and this is too small to be measured. However, violations of the equivalence principle could 
in principle lead to observable effects both for the expansion of the universe and dark matter. These would probably 
be suppressed by a small factor which is the extent to which the equivalence principle is known to be valid. 
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B. Lense-Thirring 



The Lense-Thirring effect is a gravitomagnetic effect due to the rotation of a source mass. It is difficult to measure 
and has been searched for in several experiments [66-70] but no undisputed, conclusive measurements exist yet to 
better than 0(1). Given the success of atom interferometers used as gyroscopes, it seems worthwhile to consider 
whether gravitomagnetic effects could be measured in an atom interferometer. To understand the effect, the metric 
outside a spinning body can be written as [50] 

ds 2 = (-1 - 20 - 2(j) 2 - 2tp) dt 2 + (1 - 2(f)) dx 2 + 2<f ■ dxdt (91) 
Using this in the geodesic equations (59) gives the coordinate acceleration as 

where v = . The terms proportional to C give the gravitomagnetic terms we are interested in. Outside a spinning, 
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The relevant acceleration of a test body caused by this effect is then 

GMR 2 flv 



a LT = v x 



( V x C) (94) 



Near the earth's surface, using the launch velocity v ~ 10~ 7 , the highest this acceleration can be is a LT <~ 10 _13 <7, 
which is certainly above the planned sensitivity of upcoming interferometers [10]. Though we have not done a full 
GR calculation as outlined above, following the usual guess the phase shift is kaurT 2 <~ 10~ 4 rad. Unfortunately, this 
phase shift scales very similarly with control parameters as the phase shift due to the Coriolis effect (assuming an 
earth-bound atom interferometer). The Coriolis effect gives a phase shift 

k-a c T 2 ~k- (nxvi^T 2 . (95) 

On the earth this acceleration is ~ lQ~ A g. Then, by (94), the Lense-Thirring effect is roughly a factor of <j> smaller 
than Coriolis in magnitude: a LT <~ §a c <~ 10~ 9 ac- Further, they scale the same way with the control parameters 
k, vl, and T, although their directions and dependencies on the directions of the vectors involved, vl, x, and Q, 
are different. This means the Lense-Thirring effect cannot be directly measured beneath the much larger Coriolis 
background. However, the Coriolis effect is a kincmatical effect and is thus qualitatively different from the dynamical 
Lense-Thirring effect which depends on the rotation of the source mass itself and not just of the laboratory in which 
the experiment is being performed. For one thing, it means that the Lense-Thirring effect depends on the distance 
to the source mass and not just the angular velocity. Our idea to measure the Lense-Thirring effect in an atom 
interferometer exploits this difference to isolate the effect from the much larger Coriolis background. 

To subtract off the Coriolis background, we can use a differential measurement between two simultaneous inter- 
ferometers that measure a different Lense-Thirring acceleration but the same Coriolis acceleration. Our idea for 
accomplishing this is to have the two interferometers differ in only one control parameter, their height above the 
earth's surface. This gives the same Coriolis force, up to the level that the other control parameters can be made 
equal between the two interferometers (note that time variations are not relevant here so long as the two interferome- 
ters remain identical in everything but height). However there are constraints on how identical the two interferometers 
can be. One of the most important is the need to make the launch velocities equal to high precision, since the Coriolis 
force scales directly with vl- One idea for doing this is a common launch of a single cloud that is subsequently split 
into two at differing heights. This limits how far apart the interferometers can be. We will take a height difference 
of 1 m as an optimistic but not unreasonable guess. Then the size the Lense-Thirring effect that can be measured 
is reduced by a gradient factor to olt^^ ~ 10~ 19 . This is roughly a factor of 10 4 below the initial sensitivity we 
are considering. As mentioned in Section II, there are many possibilities for improving this sensitivity by orders of 
magnitude. However this number seemed challenging enough that we have not pursued this idea further (we are 
explicitly not considering tying an atom interferometer to a telescope in a Gravity Probe B (GPB)-like configuration, 
since that has already been proposed [70] ) . It is also possible that there are better ideas for isolating Lense-Thirring 
from the Coriolis backgrounds that would allow a measurement, since the Lense-Thirring effect is naturally quite large 
compared to the sensitivities of upcoming atom interferometers. This is left to future work. 



30 



C. Preferred Frame 

There are many possible modified theories of gravity beyond those parametrized by the PPN parameters (3 and 7 
in metric (72). As one example, the PPN formalism includes eight other parameters that parametrize the possible 
metric modifications of general relativity [11]. There are many other, non- metric, theories as well. The full PPN 
metric [11] would still fall under the calculation method outlined above. 

We have not performed the full calculation to see what other effects from the PPN metric would be present in an 
atom interferometer, but we can easily guess one. Often, the new effects introduced by the full PPN metric can usefully 
be thought of as preferred frame or location effects, usually called violations of Lorentz invariance. Such effects have 
been pursued before as modulations of the local acceleration in an accelerometer on the earth's surface with periods of a 
day, a year, and so on [9, 11, 71]. While the atom interferometer can yield an impressive increase in sensitivity over the 
accelerometers used to do the previous searches (and has already led to improved limits [9]), the previous experiments 
were ultimately limited by geophysical uncertainties. Thus, it is not clear that an increase in accelerometer sensitivity 
would lead to an improved ability to search for such Lorentz violating effects without an equally improved geophysical 
model. Of course, the atom interferometer is much more than just an accelerometer, as attested to by the many 
control parameters and measurement strategies employing simultaneous differential measurements outlined above. 
Further, there may be signals from these PPN and Lorentz violating effects that are more than just modulations of 
the local acceleration. We cannot exclude the possibility that there are novel search techniques that would allow atom 
interferometers to provide stringent tests of Lorentz violating theories, but we leave such considerations to future 
work. 



VII. SUMMARY AND COMPARISON 

Relativistic effects in intcrfcromctry have been discussed before in several contexts [21-25]. None of these discussed 
specific, viable experiments for the post-Newtonian relativistic effects we have considered. Mostly, they focus on 
calculation methods instead of specific laboratory experiments. These methods are not applicable to our experimental 
setup. 

In [21], general discussion and motivation was given for considering the effects of general relativity on devices such 
as atom or neutron interferometers. Relativistic calculations for a certain type of interferometer were given in [22, 23] 
and were used to give a rough estimate of the phase shift that might arise in a setup similar to neutron interferometry. 
While it is possible that a similar estimate could be given for the atom interferometer, it would miss most of the 
important effects since their analysis does not take into account the laser pulses which actually form a light-pulse 
atom interferometer like the ones we consider here. None of these calculations can be applied directly to the atom 
interferometer for several reasons, including the lack of description of the laser pulses and the difficulty of solving 
the necessary equations for the full atom interferometer sequence in a general metric background. Essentially, these 
calculations only have what we call propagation phase. Further, this phase is not calculated along the trajectories 
that are relevant for light-pulse atom interferometry. The measurements suggested in these papers are also quite 
different from the our proposals. 

Atom interferometry is considered in [24] as a way to measure space-time curvature, which is essentially the leading 
order effect of (Newtonian) gravity. We are interested in measuring the post-Newtonian corrections. We are thus 
led to consider a specific setup in which the effect of the laser pulses and the platform on which the laser rests is 
crucial and must be taken into account. The formalism given in [25] does not give the final phase shift for the general 
relativistic effects we are interested in, though some of the effects are mentioned. We cannot use this calculation 
because the equations become too difficult to solve when all the post-Newtonian terms are kept to the order which 
is necessary given the precision of the experiment. Further, the effect of gravity on the laser pulses that serve as our 
beamsplitters and mirrors is not taken into account, and thus a relativistic prescription for the atom-light interaction 
is not included. Thus, in order to have a fully relativistic, coordinate-invariant calculation, we use our semi-classical 
method for calculating the phase shift in the interferometer. Additionally, in order to simplify the calculations, many of 
these previous papers worked in the linearized gravity approximation where the metric is expanded as g^ v w r/^ + h^. 
This cannot yield a correct result for the non-linear effects in general relativity. 

We build and expand upon this previous work by considering a specific experimental scenario for a light pulse atom 
interferometer, in which the GR effects can be calculated in a fully relativistic framework. This includes important 
effects such as a relativistic treatment of the laser pulses forming the beamsplitters and mirrors which accounts for 
the influence of gravity on the propagation of the light, as well as changes to the phase shift formulas. This framework 
is able to go beyond linearized gravity to reveal the effect of the non-linearity of gravity on an atomic interferometer. 
This requires calculating phase shifts for terms in the Hamiltonian higher than quadratic order. We also consider 
several specific experimental strategies for testing general relativity in an atom interferometer. 
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APPENDIX A: EARTH FIELD CALCULATION RESULTS 

The coordinates for point C using metric (72) as discussed in Section IV A are 
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